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"SPHERICAL WAVES IN A PLASTIC LOCKING reer 


Mario Salvadori, 1 F. “ASCE, Richard d Skalak, ASCE, and 


oo _ Asymptotic. short-time and leng- -time solutions are derived in n closed form 
for the spherical wave propagation in an infinite, plastic locking medium due | 
tovarious types of time-varying pressures applied t to the surface of a spherical | 
7 cavity. Among these are: : (a) a a constant pressure; (b) a center of dilatation; 7 
an adiabatic gas expansion; (d) an energy due to an exploding 


INTRODUCTION 


eo 
: The propagation of plane waves in various types: of locking media has been 


considered by several authors.4,5,6 Spherical waves in plastic locking media 


; Note.—Discussion open until July 1, 1961. To extend the closing date one month, a 
_ written request must be filed with the’ Executive Secretary, ASCE This paper is part 
of the copyrighted Journal of the oe Ey Mechanics Division, Proceedings of the 
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were investigated by A. _S. Kompaneets? and more recently (1960) by N. V. 
_ Zvolinskii who derived the wave equations for an elasto- -plastic medium of lim 
a Since, in general, closed-form solutions of F this type of problem cannot be 
obtained, asymptotic closed- form solutions are given in this paper for an 
finite ideal locking medium. An equation of motion is first derived for large — 
strains and displacements when a variable pressure is applied to the surface 
of a spherical cavity. It is then shown that its early time asymptotic ‘solution © 
_ is the rigorous solution of the same problem when small strains and displace-_ 
ments are considered. The long time asymptotic solution of the same equation 
is then shown to be the rigorous solution for a point source pressure when —_— 
strains and displacements are again considered. 
_ The behavior of a three-dimensional plastic locking medium with | 


under c compressive stresses is characterized as follows: 
1. At first the material of “initial density p, does not offer resistance to | 


compression, , but, once a critical v value of the dilatation is reached, such that 
the density reaches the value Py a ‘constant value of the dilatation 


wall 


1 - Po/py < <1. 


is maintained from then on (Fig. 1). 
we During the ‘Subsequent 1 motion with a constant it density p; the material 
satisfies the three- dimensional yield condition derivable from Coulomb’s law 


of 


7 San stresses, the constant c is proportional to the cohesion, andkisa ne 
tion of the coefficient of internal friction, 
The p previously described behavior is typical of certain ‘cohesive | granular 
soils under high compressive stresses. Since under such stresses the effect — 
“ofc cohesion becomes negligible as as compared to t that of friction, Eq. 2 becomes” 


as “Shock Waves in Plastic Compacting Media,” by A.S. Kompaneets, Proceedings, 
Academy of Science, USSR, Vol, 106, No.1,1956, pp. 49-52. = 
an. “On the Emission of an Elastic Wave from a Spherical Explosion in the ee.” ai 
_ by N. V. Zvolinskii, P. M. M., Vol, 24, No, 1, 1960, pp, 126-133, : : 
“On Coulomb’s Law of Failure in Soils,” by R. Shield, Journal "Mech. Physical 
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No tation. .—The letter symbols adopted for use inthis paper are defined where 

_ they first appear, in the illustrations or in the text, and are a: alpha- 

ts, for convenience of reference, in the Appendix. are 


Consider a spherical cavity of in an ideal 
«eee From t = 0 on, the cavity boundary is subjected to a uniform pressure 
w of intensity p, which is a function of the expanding cavity radius Yo. Letr be 
_ the Lagrangian coordinate of a particle, and w its radial displacement at atime — B, 
ty when the medium has been compacted to a density , py up to adistance R from © 


‘the cavity center, while beyond R, the medium is at rest with adensity p, (Fig. | 


Under t these assumptions conservation of mass requires that 


For << Be: and at an early | time when the change in the cavity radius is negli- | 
gible, so that ro ~ rj, the ratio x approaches the value : 


— 
— | 
q 
a 
« 
| 
i, in particular, that at the cavity bou ee 
in particular, that at the cavity boun 


As 
SPHERICAL WAVES 


At late times when R is large enough for (1 - to o be much less 
the ‘approaches the constant value — 


= 


Denoting time cote, the particle at r = becomes, 


and at the wave front ‘ei = R) 


ticle acceleration, by: 8, 


py 


in which w is defined by Eq. 10, ‘and its sparticular | eoletion satisfying the bound- 


ape 
-1 
|e 
o} 
- "9 @ 


and, for k = 1/2, - 


1/2)... (13b 


7) 
| 
= ‘The par 
ae ‘The equation of motion for an elementary volume of soul in spherical coor- a 
dinates (Fig. 3) is given by 
— 
‘a 
— 


February, 196 1961 
Conservation of momentum across the wi wave front (r= =R) requires, 


tion on governing t the motion a the wave front ‘becomes ieetti, 


a 
“dR ale py E(1- 


which the variable coefficients a and Bare tor by 


1-&) 


ay (xt 


(k 


= 1/2) (16d) 


Time cavity, small strain, small displace- 
ont solution) .—Substituting x = x, from Eq. 6 into Eq. 16a and b, one obtains 
the differential equation valid for small displacements. A further a 
s obtained by considering the strains to be small, that is, neglecting terms in 
# and, hence, the third term at the right hand ‘member of 8 in Eqs. 16a or 16b. y 
‘Since Eq. 6 is valid when the expansion of the cavity is negligible, and the pres- 7 
sure p is considered a function | of the cavity radius To. the pressure i is to be 
assumed constant 


an explosion!0 of. yield \ W and gas constant y 


_ “A Method of Concealing Underground Nuclear Explosions,” by A. L. oy 


LeLevier, E. A. Martinellt, and W.G, McMillan, , R-348, The RAND Comp. 30, 


— 
— 
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thus” becomes 


Long Time Asymptotic Solution (Point Source finite lowe displace- 
solution) .—Assuming for x the value £~+/° given by Eq. 7, the coefficients 
; a and £ of Eq. 15 become constant, and depend on — andk only. This equation, | 
_ rigorously valid for rj = 0, coincides withthe equation derived by Rompanseta. 7 
a The following particular solutions of this equation are of interest: ae 
a Center of Dilatation. —The magnitude ¥ of aad total force exerted on ‘the 


corresponding solution of Eq. 15 Sis 


If at a time t=t when ‘Ry, Po becomes zero, the wave front 
witha 


£2/3(g- 2a@)(1- &) \™1/ 
Since B/a 2 tow all values of & and k, R? approaches zero asymptotically 
and the wave front never stops. This result stems from the neglection of the 7 
cohesion c in Eq. 2: a finite value of c c gives a finite time at which the motion. 


in which y y and K are constant characteristics of the | gas s and of the —_ 


| 
| 
3 
sure p 4 
deal gas the 
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<3y there are no of the p point source this 
case also, the introduction of a finite value for the cohesion introduces a finite a 
_ Exploding Point Mass.—Consider a point mass M located in the medium at ; 


= 0, which at at dae =0 0 explodes 80 that its particles acquire a radial velocity Fo. 


ssuming that these particles remain on the boundary surface of an expanding 
spherical cavity y of radius Tor the t —_ kinetic energy at a time t is 


— 
| 
| 
dR 
aR 
“ate -¢ - 1) R R (8B) 


Inthe particular case of of k = 1 no other energy can occur within 


"Hence, | Eqs. (26 and 27 
“m/ dR 


pore 


solution of 1 Re 29 with R(0) 3 and ‘hence an energy input > M Rete 


“in in which < 1 for 0 
Intermediate Range Solutions .—The variable coefficients a and 8 of Eq. 15 
"approach: their asymptotic values for become practically constant 


8 versus Xg = R/rj is given in Figs. 5(a)_ ‘in 5(b) for k =1, “1/2, and 0 when 
a The solution of mq. 15 in the intermediate range between the two asymptotic | 
short-time and long-time solutions may be obtained by numerical means start- 
“ing at the limit of the range of validity for either solution and integrating for-— 
_ The constant K, needed to continue the intermediate time solution for an ex- 


-plosion is obtained by Eq. 23. 1 the range of of validity of the 


asymptotic solutions. 


Ss 


The rate of energy loss’ 
— 
— 
4 
— 
1.02 1.06 | 1.08 | | 1.25 10 | | 
7 | 1,02 | 1,04 | 1.06 | 1.08 | 1,25 | 1.50 
— | 1,00 | 1.02 | 1,04 | 1.05 | 1.07 1.21 | 1.40 | 2.14] 2.15 | 215 
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The ‘The authors are indebted to Hans H. - Bleich for the solution of the. 


PPENDIX. - —NOTATION ON 
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function of the ‘coefficient of internal friction 


= coordinate of a Particle 
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= radial displacement of a a particle 


— 
= stress 


= gas. constant _ 
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SYNOPSIS 


modes as generalized this purpose, it is frequently 
ficient to use the first few frequencies and normal modes of the building. The — 
Standard matrix iteration method is not suitable for this ‘purpose it 
yields the highest ‘mode first. A iteration method using the 


“ness matrix of ; any y shear building is s given by a and 
directly from the givenphysical data of the building. A numerical example il- 


. Ae 
BUILDING 
~The e numerical analysis of the response of an elastic rigid 
building to earthquake or blast loads presents difficulties due to the great com- 
“4 plexity of the building itself. This may | be considerably simplified by theintro- — 
duction of “shear first by P. White, 2M. 7 


2 -Note.—Discussion open ur until July a 1961, To extend the ‘closing date one month, a 
written request must be filed with the Executive Secretary, ASCE, This paper is part 


American Society of Civil Engineers, Vol. 87, No. EM 1, February, 1961. 

Assoc. Prof,, City College of New York, New York, N.Y. 
Discussion by P, White of “Natural Periods of Uniform Beam, » by 
Jacobsen, Transactions, ASCE, , Vol. 104, 1939, p. 43 


| 
The dynamic response of tall Caused by earthquake or blast load- 
| 
re 


we mean the ng simplifica-— 
1. All flo floors are assumed to be rigid. As a direct consequence of this, all 


= All columns are assumed to be inextensible, _ 
4 


3 The mass of the building is lumped at floor levels. This reduces an N- 


story building to | N lumped mass, degrees- -of- freedom problem, 


(ae.2 The effect wa direct stresses in columns on their flexural rigidities is 
eglected. 


Generally sp speaking, the first step in a dynamic- response analysis is the 
~ calculationof undamped natural frequencies andnormal modes of the buildings, 
This paper deals with a method of computing normal modes and natural fre-— 
quencies of undamped shear buildings. There are many methods available | for 
the determination of natural frequencies and normal modes, including the well- 
known Stodola-Vianello method,’ Holzer method, 4 and matrix-iteration meth-— 
od, 9,6,7 The matrix iteration method is convenient in some instances, partic- 
ularly if the work is carriedout on electronic computers, The drawback of the = 
standard matrix- iteration method is that we obtain successively the et, 4 


and lowest mode is obtained last of all, and is the least accurate of all mal 

_ Whereas in some engineering analysis, the lower frequencies areusually more 
important. Adequate results can sometimes be obtained without using all of the 7 
_ higher frequencies and modes. We show here that by using the flexibility ma- p : 
- trix of the shear building, we e may obtainthe lowest frequency and nae mode 


FLEXIBILITY ME METHOD 


inertia of girders is assumed to be 


) ey The spring constant, kj, is the shear force necessary to produce a unit rel 7 
ative horizontal displacement between the ends of all the columns of the ith 


story. Because the rotations of the ends of these columns aren not pean, 
due to the infinite rigidity of the girders, 1 we 


have 


ie, 3 “Earthquake Stresses in Shear Buildings,” by M, G. Salvadori, Transactions, ASCE, 
4 Discussion by J. A. Blume of “Earthquake Stresses in Shear Buildings,” by G. 


Salvadori, Transactions, ASCE, Vol. 119, 1954, p.171, 
:_ “Dynamics of Framed Structures, * by CG. L. - Rogers, John Wiley and Sons, Inc., New 


6 “Engineering Analysis, by S.H, Crandall, McGraw-Hill Book Co., Inc, , New York, 
i. “Elementary Matrices,” by B, A. Frazer, ua J. Duncen, and M. Collar, Cam- 


i 
4 
i; 
ee 
— 
y 
An N-story tall building is shown schematically in Fig. 1. In accordance 
with theshear building assumption, all masses of columnsare either neglected 
: 
: 
| 


SHEAR BUILDINGS 
- 4 in which E is the modulus of elasticity, in pounds per square inch; =F denotes — 
the sum of moments of inertia of all columns of the ith story, in inches to the 
fourth; a and 1; refers to the length of columns of the ith story in inches, 


1 of motion for ice mass on the sealord is 
rare 


which m, is the mass at the floor seconds squared per inch; 
and xj denotes the horizontal displacement of the ith floor, 


It is noted that « each of contains three terms on side, ex- 


cept the first and the last equation, which contain two terms. _ Eqs. 2 are used © 


— FIG, 1.—N-STORY SHEAR BUILDING 
— 


February, 1961 
by ‘most ast weitere” 9as the equation n of motion in dynamic of buildings. 
- It should be pointed out that this equation is a direct consequence of the shear 
i“ _ building assumption (infinite rigidity of floor girders), and is not valid for elas- — 
” tic rigid frame buildings, in general, If the moments of inertiaof floor girders — 
oe are not infinite, then each of Eqs. 2 will have N lefthand-side terms, Further- _ 
_more, the coefficients of are no ‘but m must be obtained 


Substituting into Eq. 2, we obtain a wee of N N simultaneous hoi ogeneous equa- 


coe eee eee ees 


It is cor conventent t to write as one matrix equation as 


[k] [5] = w*[M] [5]. 


eee 


are matrices defined as 


ae “Earthquake Forces on Systems with Several Degrees of Freedom, * by B. Hojlund ¥ 
Rasmussen, Bulletin of the Seismological Society of America, Vol. 45, October, 
9 “Response of Tall Buildings to Random Earthquakes,” by A. C. Eringen, Procee - 
ings of the 3rd United States Natl. Congress of Applied Mechanics, 1958, p. an 
10 “Stiffness Method of Rigid Frame Analysis,” by M. L. Pei, Proceedings of f the all 7 
ASCE Conference on Electronic Computation, September, 1960. 


ow 11 «The Dynamic Response of Tall Structures to Lateral Loads,” by L. ee, Pro- 


| 
- ky 5, + (Ky + kg) 59 - kg 03 = — 

| 
in which [ 

7 
g ima 
a 


Mass matrix x [Me 


Premultipiying Eq. 5 by the inverse of [ K], yields 


6 


which, for simplicity, we [F] =[ K] - and [H] =[F]({M]. If matrix 

_ iteration procedure is carried out with Eq. 9, we shall obtain the fundamental a 

- mode and natural frequency of the shear building. Usually, matrix —— 
is a lengthy process, and it is not too practical to calculate Eq.9 for a sore 
building. However, , for shear the exact of K is given by the 


formula 


y=) 


which. fj is the reciprocal of Or, written ‘out in full, 


fy 


— > — 
- 
— 
: i 
— | 4 


iL __ The proof of Eq. 10 is fairly simple, ‘since the elements in the flexibility 
: _ matrix are the displacement influence coefficients of the shear building. That _ 
is, fjj is the displacement of the ne ith | floor due to a unit load applied at the jth 
7 7 floor level, With the unit load applied at the first floor of the peep, shown | 
( in Fig. 2, the displacements of all floors arethe same and equal to Ly 3/12 E a 
Hence, every element i inthe first column of the flexibility matrix is i es 


FIG. 2. —DISPLACEMENTS OF 
TO A UNIT LOAD AT 
TH E FIRST FLOOR 
LEVEL 
"Next, with the unit load applied at the second floor on. the | displacement of 
Bars Sret floor i is equal to f;; and the displacements of all upper floors are _ . 
Thus, the second column of the flexibility matrix is, 


_ The remaining columns of the flexibility matrix canbe dovtwetin the same way. | 
Briefly, the iteration procedure is as follows: 


- Compute the flexibility 1 matrix x directly from given data by means 0 of Eq. 


2. Assume arbitrarily a column vector as the initial 


is the al mode, 


_ in which '™ and Or+1 are the corresponding values of 5 of the same » floor in 
; i the rthand (r+1)th iteration, Usually, we may use the 6 having the largest ab- 
solute magnitude for this calculation. Some idea of the accuracy y of the result 
7 : may be obtained by means of the enclosure e theorem, 6 which gives the upper a 
and limit of the true frequency as, 


a better estimate of the may be by taking 


the ratio of the weighted average eof and — 


4 

ee ee 


i= 


‘ 


First Mode. —Although the method is particularly for tall 
we shall for 5 sake use a tl three- -story building analyzed by Rogers” 


ALYSIS DATA 


Mass, in Ib~sec” per in, 


- 
j 
= 
— 
& 
4 
= 
- 
; 
Spring constant k, in kips per inch 
— 


‘February, 


as an illustration, The building is shown in Fig. 3, and its masses ond owing 


nstants are given in Table 1, The mass matrix is a diagonal 1 matrix, al 


‘fito 3 36 36 


went 
in,per | 


For? 


oa Since we are concerned with the “ilies magnitudes of the displacements 


of the normal n mode, we we may discard the numerical factor from H to simpify 
calculations. To begin the calculations, we assume 


oO 


36 


43 8435 
= "660, 780 


1,258, 180] = 660,780. 


2.48027] 


— 
cycle.— Applying Eg fio} 
lao 


LDINGS 7 
‘This may may = be taken as the first normal mode of the shear building, The natural 
frequency is obtained from ~ 


5, o/ =4 = 10, 56 radians per sec 

w 


= 11,06 radians per 
@g= 11, 12 radians per s¢ sec 
11,13 radians per sec 
Alternatively, the fundamental frequency may be obtained by Eq. 1 13, Using, for 


ome, the results of the first and second cycle, we have > 


4 


= 11. 13 radians per 
(3420? + 6500 + 84352) (15/360,000)? 


. — Having obtained the first mode, we may proceed to calcu- 


late the next higher mode, The general procedure is well known$, 7, 12 12 and will 7 
‘not be discussed here. In this Paaee, one may ¢ calculate a as many, or as | few, 


BUI 


A Method for the “Solution of Matr ces,” P 
Vol. 7, No. 17, 1934, pp. 865- -909. 
2. «Numerical Analysis o of Earthquake Response of aTall Building,” by 
oq ‘Tung and N. M. Newmark, Bulletin of the Seismological Society of America, 
Vol. October, 1955, p. 269. 
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«Notes on the Determination of Higher Modes of Vibration by the Stodola o 
Iteration Method,” by H. E, Fettis, Journal of Aeronautical Science, Vol. 26, No. 
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INTRODUCTION 


_ One of the major problems facing structuraland mechanical engineers is the | 
design of structures which are capable - withstanding Leap pressure loadings 
of a transient or time dependent nature. For example, in the rapidly growing 
field of protective construction, the design requirement that structures should - 7” 
withstand high intensity overpressures produced by nuclear blasts, may lead for 
economic cand other ‘practical r reasons to the anager of abandoning the conven- a 


working stress concept). Instead, plastic flow at various locations in the struc- 
ture is permitted, but the total deflection must be kept within acceptable limits. _ 
An elastic- -plastic analysis of such structures under transient pressure 


required to determine two major design 


The maximum permanent deflections in structure due to. its 


a Note.—Discussion open until July 1, 1961. To extend the closing date one month, a fl 
| written request must be filed with the Executive Secretary, ASCE. This paper is part rs : 
of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the 


American of Civil Vol. 87, No. EM 1961, 


ee 1 Assoc., , Paul Weidlinger, Cons, . Engr., New York, and Adjunct Assoc. Prof, of Civ, 
2 Dir., Inst. of Flight Structures, Columbia Uni Univ., 


Cons. >. Engr. New ‘York, 


At 
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ious could be formulated by writing the differential equations of motion - 7 
a three-dimensional elastic-plastic continuum. The complexity of the resulting 
_ equations for the complicated geometric shapes involved would obviously make 
their solution a major problem; even ‘in the simple elastic problem, the treat- 7 
* ment as a continuum is only manageable in a few simple cases. As a result, 
the beam concept has been introduced, replacing the structure by a center line 
with: stiffness properties (that is, a linear relation between the bending moment 
and curvature). This approach gives good results in the elastic caseas long as 
F a length of the structural member is s sufficiently greater th than the transverse 
_ dimensions | of the member. 


2 
uJ 
= 
2) 
= 


BENDING 


— 


similar approach is used in the of elastic-| and 


_ structures. A relation, which can be verified experimentally, between the bend- — 
- ing moment Mand the curvature of the deflected center line 1/R is the basis for : 
_ all current elastic-plastic analytical methods. This relation has the character _ 

shown in Fig. 1 and consists first of a linear portion followed by a transition — 

Paiow a nearly horizontal line. For the sake of simplicity the actual bending mo- 

1 ment-curvature relation is replaced by an idealized curve consisting of the two Po: 

_ straight lines shown in Fig. 1. This implies that the curvature is restricted to - 
those values for which the strain hardening range is not entered. It also means 

that some detail near the knee of the M - 1/R diagram is lost, however, for 
members of structural steel, the ideal diagram is quite close to the actual one 


and the differences | in the result must therefore be small. It should also be noted — E 7 


‘The dynamic and the load transmission through the struct 4 
me oonto i ing foundation may in many cases be a controlling factor in a 
ae ie elastic beam or frame structure which is su 4 ra 
| 
: 

= 
aa 
yf 
— 
™ 


may be useful and appropriate for special 
In structures es subjected to loads en produce large permanent deforma- 


perigee to the energy absorbed in n the plastic stage and consequently _ 
the elastic response of the structure might be entirely neglected, leading to a 


2 a . In structures for which the location of yield hinges can be clearly pre- 
dictedfor given pressure loadings, an with with plastic hinge locations which 


assumed a priori | can be performed. ad 


Botha assumptions result inattractive analytical the 
_ propriateness of these methods, with the exception of certain trivial cases, re-— 
mains ‘somewhat in doubt unless they can be compared with more general pro- — 
_ cedures which do not put undue restrictions onthe location of the plastic hinges — qa 
_ or on the duration of the successive plastic and elastic stages at the points of : 
formation on the structure. In recent years, several less restrictive theo-_ 
ries of dynamic elastic-plastic analysis have been formulated, but practical ap- 
_ plications have been limited to = structures, due to both mathematical and 
A theory was presented* 1959 which permits a general analysis in 
i.) ficient detail for most practical problems, and its application by finite difference — 
methods has been found to be ideally suited for rhigh Speed computing 
7 


1. To present in detail the finite difference method for the 5 Niall -plastic | 


analysis of structures, 


— clarify the appropriateness and range of validity of the simpler theo- 


a results of the finite difference method. 


is used, in which po the static capacity load which will “produce 
y ‘. yielding in the structure under consideration; n Po. is the initial value of the dy- 
‘namic load at the time t = 0; and t, is the decay constant. While the finite dif- 


ference method ¢ can be generalized for to complex- -frame aye struc- ay 


- ries by comparing the results obtained from these theories with the more . pre- - 


= is given in this paper and examples will be presented for 
simply supported beam only. The range of validity of the simplified theories — 
_will be generally defined by the intensity | of the initial value of the load interms | 


: 4 of the amplification factor n, W while | uSing the value of the maximum permanent > 
defile 


deflection of the beam a asa ‘comparison criterion, 


“Dynamic Design in the Plastic Range,” by H. Bleich and M. Baron, ‘presented 
at the 1959 ASCE Meeting, Mechanics Div., Los Angeles, Calif. 
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ay Notation. —The letter symbols adopted for use in this paper are defined where — a 

fe 


they first appear, the illustrations or in the text, and are 


FINITE METHOD FOR THE RESPONSE OF ELASTIC- 

_ PLASTIC STRUCTURES TO TRANSIENT — ae 


nto a series of normal mnenen of vibration was presented by Hans H. Bleich, :. a 
F. ASCE. In this analysis, the original modes" of t the elastic ‘system were vutil- 


A method of analysis involving an expansion yn of the response structure 
i 


Frown of the formation of multiple plastic hinges (hinges at various i 
, on the structure) and for the re-elasticizing and re-plasticizing of these points. 
Utilization of some of the general principles presented by Bleich led to. the 
_ development of the method of analysis of this section. * The method | which em- a 
_ ploys a finite difference approach is presented in detail for application to stat- 
ically determinate structures and is particularly well suited for use on high | 
- speed computers such as the IBM 704 which was used in the solution of the pres- 
problem and other problems.6 = 
7 ar _ At this point, one may consider the type of problem which can be treated by 
anormal mode, or by a finite difference approach in general, with or without 
—_— . Instructural problems, the designer is usually satisfied with aknow- 4 
ledge of data on a beamat points where the Spacing is of the order of the beam 
depth or perhaps at the most half of the beam depth, and both modal and finite. 
: difference analyses are suitable for such problems. . Mode methods are useful 
because computations | be restricted toa relatively small number 
_ modes and, similarly, finite difference methods are permissible since the com- 
_ putations canbe restricted toa relatively small number of mass points, leading 
to sufficiently accurate, physically, meaningful results. The elements of the 


_ structure will be assumed to be ideally elastic- plastic. ‘They will exhibit the ’ 
- conventional elastic- -plastic behavior of bending, that is to act elastically up a ; 


- reached (Fig. 1). When the capacity moment Mc, is maintained ata particular 
section, curvatures of arbitrarily large magnitude can occur there, and plastic- 
"Since the curvature at a a section theoretically approaches infinity as the mo- = 
ment at the section approaches M,, one may postulate that the attainment of the 
fully- plastic moment ata section allows finite changes of slope to occur a at that 
; Ss ection while the moment remains constant and equal to Me. This corresponds ~ 
. £ to inserting a physical hinge which can only transmita constant moment M, at 
§ that section, hence, the concept of a plastic hinge. The capacity m moment Me of 
a member may be selected to allow fordirect stresses whichoccur - simultane- 
ously with the moments, if these stresses are large. 
_ Equations of Motion-Elastic Stages.—Consider a simply supported beam of 
ui length L (Fig. 2) which can be replaced by a continuous chain of equally or un- a 
_ equally spaced lumped rigid masses connected by inextensible bars with mo-- 


_ ‘ment springs at each _ joint. Let E refer to the modulus of elasticity | of beam 


3s “Response of Elastoplastic Structures to Transient Loads,” by H. H. Bleich, -Pro- 
ceedings, New York Academy of Sciences, November, 1955. 
i. “On a Method of Amplification of the Dynamic Strength of Structures,” by P. Weid- __ 
‘linger, M. L. Baron, and R, Skalak, -— Air Force Special Weapons” Center, Report No. 

-TR- 10, October, 1958, 
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a FINITE METHOD 
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Me 


2.—BEAM MODEL 


‘FIG, 3, —NOMENCLATURE FOR MODEL 
ite 
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4,—-DETERMINATION OF 6, 


_ The equations of motion of the masses during the —_ ic | 
“rived by utilization of Hamilton’ principle and the calculus of variations. 
_ Using the nomenclature of Fig. , 3, in which my is the mass of the structure 
er unit of length in the span hn, the lumped mass at | the point | n, tm, ~ defined 


m 


and the total energy ry of all the mass 188 points of of ‘the is ‘obtained by 


n n+1n+1| , 2 


7 


is thet moment required to | produce a unit- -angle rotation at the joint The 
is obtained from Fig 


=- 
n 


h 
‘Using 7 and 8, ‘Eqs. 5 and ‘beco 


me; 


n-1 


a 6s material, I be the moment of inertia of cross section of the beam (constant ae “4 

i all sections), m denote the mass of the beam per unit length of span, andy re- et 
on 
ae 
tm 
The kine W = 7 = 

The potential energy of the structure, that is, the strain ene 

ts 
| 
— 
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FINITE DIFFERENCE METHOD 
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EW 


we obtain, by ‘means of the calculus. variations, ‘equation of motion for the 


free vibrations of the th concentrated 


h 
ast 


the structure is subjected to loadings Qn ( (t) on each span hh wt t 


"3 
a 


> 
= 
— 
— 
4 
Ve 


al (Note that q, dn is the load p per | ‘unit of in ‘the In the 


points are equally spaced on either side of my» that is h, + k= =h (k = 0, 1, = 
the equation of motion for the concentrated mass m, becomes 
m +m. 1) 


cat 
(19) 


oe 


and using Eq. 1 


Referring to the beam in Fig. 2 | (b) wnere itis is wn into ten equal intervals 
th h = L/10, the structure is analyzed asa ‘continuous chain of eleven 
_ lumped mass points located at the positions 0, 1, 2, 3, 4, and 5. As the loading © 
onthe structure will be symmetrical with respect tothe center line, in the pres-_ 
ent example, only | the equations of motionfor six of the masses will be required. 


‘ Eqs. 18 and 21 may be used for the concentrated masses atthe points 2, 3 


of leng 


— 
7 
— 
a 
— 
4 
pee or a constant moment of inertia, the equations f [| become, re-_ eae 
| 
| 


a enna: of the equations \ of ‘motion in the elastic range for the model of 
Fig. 2 (b) is as follows 
2 M - 2 M ] = 2a 


<> 


= Replacing the function Yn (t) by the first term of its finite differ difference ex cere 
pansion in 1 central differences, palin ae 


(t+k)-2y (+y, 


in which k is the interval of the time steps, recurrence equations 3 for the de- ¥ 
; flection y Yn (t+k) are obtained in terms of the moments, loads and displacements = 


‘ (250) 


& 


— 
in which the moment-displacement equatioms are piven 
— 
| 5 
4 
— 


The initial conditions of zero and zero 


a and to the starting formulas 


trols the choice the time step “k” is examined in Appendix I. 

a Plastic Stages .— —For ‘sufficiently large pressure loadings, Gn (t), plastic de- 

‘eceumionse may start at a timet = =t at any of the various points n, on the beam. 

ee deformation curve of the structure will no longer be smooth, 
4 


but will show “kinks? at ‘points: where the plastic hinges occurred at Sg? time 


The question of the stability of the of | 25 con. 


time t the moment M. computed from Eqs. 23 will be larger the 
_ pacity moment of the beam, Me. This is physically impossible and so a cor- 
- rective displacement must be introduced at each point on the beam in which M, : 
</ Mj /, Such that the condition /M,/ < Mg is enforced at every point on the 


corrective displacement due to at | the joint j only is s 
given 


in n which Yn isthe displacement at the joint n due tothe corrective displacement 7 
plasticity at the joint is a coordinate proportional the 


§ carrer n of a fictional elastic beam ifa hinge were introduced at the point x = xj 

7 and the two portions of the structure on either side of the point j were rotated 
acertain angle. This angle is selected so that is equal to minus unity 
(of length), in which case the plastic kink angle at the joint j is equal to 6; j/he : 
displacement is shown in Fig. 5 for the case of the simply supported 
“The deflection curve consists of straight | lines witha kink at the point n = jwhere 4 
the plasticity occurred. “plasticity” occurs at several points, the 
displacement at any joint n due to the plasticity at tthe s several a , is is ob- — 


— 4 
y, (0) = 0, (kK) = y, (-k)... (27) 
: 
: 
5 
; 
Po. 


expression for the: true deflection a at point is by 


in which Yn. isthe vat the joint 2, A 
q 


tion for the moment M,, becomes 
Writing Eq. 32 32 as a 
a 
ve. 


KINK ANGLE SUCH THAT 


ten without the summation 


4 
&g 

&g 
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in 
(By definition 6 Won = 1, in order that the kink angle en 1 be senttive. ) 
The plastic angle “generalized coordinate at the joint n, On (t) may be — 
ten in terms of the previous plastic angle coordinate at the time stept - -k, plus © 
the in the in the intervalt-ktot 


‘ 


in which 


Once e deflections Yn (t) have found from 25, Ea. $8 


me compute t the moment t M§ (t), and the e condition that that t /M nt Mc ‘is then used 


to evaluate the change inthe plastic angle, coordinate, A a(t - =k). The moments 
My (t) are then evaluated by means of Eq. ' 37 and are used in Eqs. 25 for the de-- 
termination of the deflectiony, (t+k), 

It remains now todiscuss the procedure for finding 46, (t- The incre- 
“ment of plastic angle may be positive, thus indicating an increase in the plastic f | 
deformation of the beam at the joint n in the time intervalt -k tot, or zero, 

thus indicating no further plasticity in this time interval and a return of the aa 
joint n tothe elastic range. It should be noted that when a point of the beam. re-_ 
enters an elastic Stage, it maintains a permanent plastic kink with an angle 

_ proportional to @n = 9p (t - k). Moreover, any or all points on the structure 

may enter and leave the plastic stageat any time with no fundamental difficulty 

in applying the equations. 

_ To find 46, (t-k), let the plastic angle e Onb be e known up to the time step 
t - k and consider the increment of plastic angle in the interval t -ktot. The 
procedure canbe used: 
; 1, Compute the deflections yp () from the equations of motion of the mass 
Using the values of y,, (t), compute the moments M,, (t ) from Eqs. con- 


— 
‘ 
7 


ppeve if the structure were to ian the deflections Yn “O plus the effect o Pe 


all previous plastic deformations up to the time step ae 7 ; 


_ 4. Compare ME (t) withthe capacity moment of the sian’ Mc (t). If /M 
_K Mg, then the beam has reentered the elastic stage a at the point n and there — ; 
is nofurther plastic flow at this joint in the interval t - k to t. If /My, (t)/> My, ; 

- i“ additional plastic flow occurs and the change in the plastic angle “coordinate oF 


‘The coordinate is then 
~ 


moments which are ‘to be used in Eqs. 25 to ) obtain the the deflections 
the beam in the time Yn (t +k), 


General Methods.—A general methcd analysis of | elastic- -plastic structures” 

in which plastic flow may develop at various points in the structure, was pre- 
Be by Bleich and Mario G. Salvadori, F. ASCE. This method uses a nor- : 

mal mode of vibration approach; one a is used during the elastic stage 
and when the first plastic hinge appears, a new set of modes is used. When- iy 

ever a new plastic hinge appears, ora significant change in the location of a 

hinge occurs, different modes must be used. This is also true when a section 
or _The is s effective in in relatively — 


‘The method is more complicated problems where many sig- 
nificant plastic hinges occur, the determination of the modes and 
changes required to go from one set to another require considerable computa- pS 
tional work. Moreover, convergence problems v with the series of modes for c- 
_ termining the flexural moment make it extremely difficult to apply the method | 


to structures which are loaded insuch a manner that they reelasticize and all ° 


_ Impulse Motion of Elasto-Plastic Beams,” by H. H. Bleich and M. G, Salv 
‘Transactions, Vol, 120, 1955, p. 499. 


— MODE APPROACHES I } OF ELASTIC-PLASTIC 
a 
| 


a of A general mo de approach which of 


For simple structures in which the location of a single significant | niastic o 


ther simplified version to be subsequently presented. 
«Simplified Mode Analysis - One Mode Elastic Stage Followed ess a Rigid — 
Plastic One Hinge Stage.—For relatively simple structures under appropri- 
ately simple time- -dependent loadings (such as the simply supported beam loaded ~ i 
by a pressure which decays exponentially in time), a simple approximation to 
_ the general-mode method previously given has been widely used, in which only a 
P a - a Single elastic mode is followed by a rigid- plastic one- hinge motion. This ap- 
proximation, which be calleda “ one of freedom” method, is 


at the center of the fundamental mode of a simply sup- 


q 


HINGE 


FIG, 6.—SIMPLIFIED MODE ANALYSIS (III-b) 


- — beam with a yield hinge at the center is used. This plastic mode cor- 
responds to a rigid body rotational motion (zero frequency) and hence tothe ri- 

The conditions on velocity and displacement at the end of the elastic stage 
serve as initial conditions for the start of the rigid-plastic stage. Due to the - 4 
fact that the motion of the simple beam is approximated in both the elastic and = Gi 
- plastic Stages by a single fundamental mode, the conservation of momentum re- = 
‘quirement necessitates a jump in the velocity of magnitude 4/7 when determining» a 
the initial velocity for the plastic stage, 

“YW Computed results obtained by applying this simplified theory tothe analysis — 
of the simply supported beam will be ome in the subsequent section of this pa-- a 
“Rigid Plastic Analysis. —The dynamic analysis of elastic- ~plastic systems 
under loadings which produce fairly large plastic deformations may, under cer- :, 
2 tain c circumstances, be made ina sufficiently accurate and quite simple manner — 


ae 

we 

| 
— 
— 
“2 
| 
= 
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by assuming the system to be composed of an ideal rigid- plastic 1 material. The 7 
"| ‘solution thus obtained neglects the elastic part of the displacement and the re- 
4 ‘sponse predicted by this type of solutionis accurate only if the plastic deforma- 
_ tions are large compared to the elastic deformations. io Lee and P. and Ss. 
4 _ Symonds, F. ASCE have presented” 8 a general rigid plastic beam solution for ‘ 
loadings, while a simplified one- rigid- plastic 


DECAYING PRESSURE—FINITE DIFFERENCE METHOD 


| 


Fundamental Period of beam 
7 Effect of exponentially decaying pressure 


elastic deflection due to Po 


Total Plastic Angle 6 x 
Radiang 


= 
. 


= 


0.0195 


(0.048 


0,034 
0,022 


{ 
= 


which also includes the e effect of shear ow was presented 9 | by Salvadori — 
Paul Weidlinger, F. ASCE. The formulas for the rigid-plastic one-hinge 

analysis of the simply supported beam under consideration are given in detail 
by Sa. Salvadori and i and Weidlinger. ’ Comparative results of the analysis and the an _ 


Large Plastic Deformations of Beams Under Transverse Impact,” by E. H. Lee 
and P, S. Symonds, Journal of Applied Mechanics, Transactions, ASME, Vol, 74,1952, 
"9 On the Dynamic Strength of Rigid-Plastic Beams Under Blast Loads,” by » M. G, 7 
and P, + Proceedings, ASCE » Vol, § 83, EM 4, October, 1957, 
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NUMERICAL EXAMPLES AND COMPARISON OF RESULTS oko 


n > fini ethod is applied to the 


simply supported beam under the in a3 


«PARUON OF 


Finite Difference Method 


Simplified Analysis (One Mode Elastic 
State Followed by a Rigid-Plastic 
One Hinge Stage) 
- 
The total plastic angle 6; is given at various points on the beam, i = 0, Es 3, 
and 4, The maximum total midspan deflection due to the dynamic load p (t) is 
given in terms of the maximum static elastic midspan deflection A, due to the - 
le’ capacity pressure Po» ‘It is of interest to note thatfor allvalues of n, the largest — 
plastic angle is produced at midspan; this angle is a than 
the plastic angle at the other points on the beam. a Sia 
_ Table 2 shows similar results for aconstant step pressure loading of ‘than 


1.00, the deflection of the beam at midspan increases indefinitely with time. _ 


nd Np, for the values n - 0.50, 0.75, and 1.00. For values of n greater than 


y > alysis considered in the two preceding sections of this paper are given inthe ~~ oe 
followin 
in which the decay constant is tp = T/2, T being the ont: 
fi : elastic beam. Numerical results for various values ; i 
Total Plastic Angle 6, x 10* in radians Ratio 
= 060 | 09 | | @© | 1,024 — 

— 
Factor ne rs + 3 — 
= — — 

| 

— 
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FINITE DIFFERENCE METHOD 


of Results.- t the ‘exponentially decaying 1 loading, tl the results 

obtained fr from the finite difference analysis are compared in Table 3 withthose 
obtained from the following simplified 

re simplified elastic- -plastic analysis consisting of a single mode elastic = 
* followed by a one hinge (at midspan) rigid plastic : stage, that is the first 7 
elastic- -plastic mode. (The simplified mode analysis presented in the preceding a. 
A single hinge (at midspan) rigid analyais. (The rigid plastic an- 
of the preceding section). 


ore LY DECAYING PRESSURE: n=2 
= -¥/ te 
q 


FINITE DI! PRERENCE 


SIMPLIFIED ANALYSIS. 


MODE ELASTIC PHASE 
FOLLOWED BY A ONE-HINGE 
RIGID PLASTIC STAGE ) 


RIGID PLASTIC 
__ONE- HINGE ANALYSIS 


‘Typical deflection-time curves for the three cases are shown in Fig. 7 for 

the case of the amplification factor 2. Fig. 8 shows a curve of Max/A 
Fa n; this curve is obtained from the finite difference analysis results. " 

Points for n = 0. 75, 2.00, 3. 00, and De 00, computed by’ the approximate —. 


—The numerical results that the authors have obtained 


from the finite difference method (Tables 1 and 2) have shown that a time de- a 


_ pendent loading, , uniform in space over the beam, produces a single well et 


q point of large yielding at midspan. At this point, a major portion of the plastic 


| 
= 
FIG, 7,—TYPICAL DISPLACEMENT-TIME CURVE g 
x 
— 


to the final deformation of the beam. “This implies that a single plastic hinge 
can be assumed to form at midspan and a simplified one-hinge analysis will 
predict | the e maximum | deflection with sufficient precision for many practical = 
_ purposes. In the particular cases that have been examined, the simplified an- q 
alyses result in deflections which are smaller than those obtained by the more “4 


7 2 _ elaborate and accurate finite difference analysis, which includes the effects of 
higher modes. It is not certain however that the simplified analysis willalways , 7m 
For more complex structures and loading situations, the simplified methods _ a 


may ~ — bad results. The authors have found that even for the relatively 


ig 


FINITE DIFFERENCE 


ELASTIC- 


on 
ONE Hing HINGE RIGID - 


4 


a +5 5.0 


nae simple case of a beam with cantilevers, the simplified methods may give re- 
= with significant errors. Generally, when a more complex distribution in 
_ ‘Space of the load is present, or where the strength of the structure is not con- 
7  stantalong its centerline, higher modes of the structure may produce significant 
contributions to the deflection. In such cases, the elastic-plastic mode analysis 
becomes enemas and the more owerful finite difference method must be 


energy is absorbed, and consequently, additional yield hinges which maydevelop i 
q 
| 


—STABILITY OF THE NUMERICAL INTEGRATION 


T 


he stability of the numerical integration in 1 time of the eq equations s of ‘motion, — 
a Eqs. 21, 22a and 22b, is considered with a view towards determining the maxi-_ 
= value of the time step k for which the solution will not divergewith anin- J 
_ Eq. 21 is conveniently written interms a the displacements alone | (Eqs. 15). . iF 


a 
these equations of motion of the beam in Fig. 2 is equivalent 
meeiwed analysis of the lumped mass beam having five symmetrical modes of 


ration. can into ‘set five of form 


is a constant, and is a linear combination of the loads dn 
The stability criterion for the solution of this system of equations, 10 using the 


“For t present problem, ‘the shortest ‘period, % of the the lumped 
masses controls. This period is approximately equal to the period of the 
_ symmetrical mode of vibration of the beam with distributed mass. - 
_ Astable solutiondoes not imply accuracy, which is also dependent | onthe ac- 


cumulation of round- off errors in the various stages of the analysis . A con- 


- siderably smaller value of k than that given by Eq. 44, bie a — be 
used in the numerical computations. — 


| 


¢ 


odulus- of elasticity of beam material 


10 “Numerical in Engineering” by M. G. Salvadori and M. L. Baron, Prentice 
He, 2nd Ed,, 1955, pp. 121-132, Problem 3.99, 


‘a . used. This method has been applied successfully by the authorsto problems of [i 
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oment 
moment 


mass per unit length sailed span 
= initial value of dynamic load at tim time t = 


= static capacity 'y load which will produce 3 sali in ae structure 
period of ‘elastic be beam 


= at joint n n due to corrective. or? 
static elastic ic deflection _ 


a generalized coordinate proportional to to the angle of ‘the plastic kink 


"deflection curve at the pc points n m a fictional elastic beam ifa hinge 
oe were introduced att the point x = x; and the two portions of the struc- 


of the point were rotated by a certain angle 
Roi ‘of vib vibration of of the n mode of the lumped mass system > 
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Axisymmetric | bending of a plate of finite supported 
_ by a viscoelastic foundation of infinite extent is examined. - Assuming friction- 
a -_ less interface conditions, a formal solution of the problem is obtained by the 
_ use of iterated integral transforms. An example involving — method of of com- 


INTRODUCTION 
> The ren a plates a a deformable foundation rests on on assumptions 7 
- concerning the behavior of the plate-foundation system. ‘These assumptions | 
4 involve: ( 1) System of equations by which the plate behavior is defined, (2) 
description | of the foundation, that is, whether a solid continuum or a set of 
——: or partially interconnected individual elements, and (3) condi- 


theories of have appeared. "Among these may be cited the pa- 
‘4 pers of Hertz for a thin plate ona Winkler foundation, D. L. Holl for a thin 
_ plate on an elastic half-space, Szabo for a thick plate on an elastic half- -space, 

_" P. M. Naghdi and J. C. Rowley for a thick plate | on a Winkler foundation. 

A short bibliography of'such works has been presented by S. Timoshenko and 
Woinowsky- Krieger.2 Avoiding a detailed review of previous work in this 


written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the 
American Society of Civil Engineers, Vol, 87, No. EM 1, February, 1961, 
Assoc, Prof, of Civ, Engrg., Univ. of California, Berkeley, Calif. 


 Note.- —Discussion open until July 1, , 1961. To extend the ) closing date one month, a 
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area, it may be 7 that the thin plate theory ( Poisson- -Kirchhoff) 1 ‘ inade- 7 
quate for treating problems involving localized surface loads. The thick plate in 
theory within the framework of three-dimensional elastic theory involves a 
cumbersome ‘mathematics, and the Winkler foundation does not permit ex- 
amination of stress conditions in the foundation itself. 
_ his paper treats the problem of axisymmetric bending of a thick plate of 
infinite extent resting on a semi-infinite solid. Iho order to retain the analyti- 

| cal advantages of a two-dimensional plate theory, | the theory of E. Reissner, 

which the effects: of transverse normal and shearing stresses, is 

adopted. The Boussinesq: solution for surface loading ofa half- space along 

- with the assumption of a frictionless interface completes the description a 

the plate-foundation system. In order to account for the possibility of visco- fay 

elastic material properties in the plate and foundation, the well-known ry value 

a spondence principle relating linear elastic and viscoelastic boundary value 

 Notation.—The letter symbols adopted. for use in this paper are defined 

3 “where they first appear, in the illustrations or in the text, and are arranged ; 

alphabetically, for convenience of reference, in 


"The equations describing Reissner’s of of on an elas- 
tic foundation of the Winkler type have been extensively presented else- : 


_ where. 4,5 Here, will be to restate the necessary in a 


id & | 


= re) 


er’ 


form: ‘suitable for application. to the problem of bending o of. 


plate resting on an elastic half- “space. . Using notation defined by Timoshenko | 


ia 3 “On the Bending « of Axially Symmetric Plates | on on Elastic Foundations,” by =, 
Naghdi and J.C, Rowley, Proceedings, First Midwestern Conf, on Solid Mechanics, Univ. 
os — “Thick Rectangular Plates on an Elastic Foundation,” "by D D, Frederick, Transac- 


“On Some Problems in of Thick Circular on an n Elastic Foundation, 
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_ VISCOELASTIC PLATES 


in which = 


=a equations w w “represents the vertical | displacement * the plate’ s middle 


at plate-foundation interface; h denotes the plate thickness: and pu tetrad 
Poisson’s ratio. Anticipation of the time dependence of the deflection for vis- 
coelastic problems accounts forthe retention of partial derivatives in the axi- 


(ne +22 te 


Mg= - w, w 
| 


stresses Or, transverse shearing stress and transverse 


normal stress o,, are obtained from My, M pe. q and and p in the usual manner 


a wo conditions associated with Eq. 1 have been examined previous-_ 


: ih ‘5 In the present problem for a plate of infinite extent, the middle sur- a Pe 


‘face deflection and appropriate combinations of its derivatives must vanish oy ws 


FOR THE ELASTIC HALF-SPACE 


complete the of the plate bending problem, an expression 
for: the surface displacement of an elastic half- -space subjected to an axisym- 
metric normal pressure is required. For the present purposes the — 
= the problem is most conveniently expressed in the form secret - Tereza- 
6 Referring to Fig. 2 the surface deflection ss 
“Hie 


6 “Fourier Transforms, by I. N, Sneddon, McGraw-Hill Book Co., Inc., New York 
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f- p(r) r J, (ar) dr. 
the Hankel transform of zero order of the surface p(r). 


stresses and displacements in the interior of the half- ~space are expressible 
in a similar form. Because attention is focused, herein, on on plate behavior, 


6 - anes 


the reader ris r referred to the work of IN. Sneddon for further details con- 
cerning the computation of stresses and displacements in the half- -space, once 


the surface ‘pressure is determined. 


_ ELASTIC PLATE ON ELASTIC 

_ Fromthe previously stated equations, it is possible to obtain axisymmetric 

: solutions of an elastic plate-foundation system. A frictionless interface con-_ 
dition is assumed, and the plate middle-surface displacement is set equal to 
- the surface — of the foundation. Operating on Eq. 1 with the Hankel 


al 
3) 
\ The 
= 
— 

— 

= 


a(a), are obtained in the same fashion. To obtain Eq. 9 have 


"assumed the. following properties of the solution: rw, rV“w, 
r ,~ tend to zero as r approaches : zero or infinity.°¢ 6a A comma followed | 


r Vw 
r 
by a subscript indicates partial differentiation \ with — to the ‘subscript. 
_ +6 h? a2 


Substituting Eq. 11 into Eq. 4 


‘The ‘Hankel integral Eq. 10) gives another 


for w(r) 


 Inasm Inasmuch a as. Eqs. 12 and 13% are fora all values r, the 


expressions mi must be be equal. 


“Inversion of Eg. 14 yields 


k 7 D 3 | 


this ec ala), the 
obtained from 4 
wall | 


At this stage it is convenient to introduce th | following notation to non-— 
-_dimensionalize the system ‘Parameters: 


Bat 15 th then assumes ‘the form | 


« 


EM1 a 
the Hankel tr - ankel 
a | a 
4 
(15) 
lied load on the plate is 


which is and thickness parameter the formal solu- 
tion for the deflection of the plate is complete. The Hankel transform of the © 


interface pressure can und by § Eq. into Eq. 11. The result 


‘Returning now to Eq. 5, a solution for the tile resultant Qy can be ‘obtained 
= Substituting Eq. 20 into Eq. 5 wii interchanging the order of — 


Toe the | bending from Eqs. and 4, is needed. ‘Differ- 


Eq. 21 and simplifying gives 


Appropriate differentiation of Eq. 18 and substitution into Ea, od gives the radial 


— 
ian 


lar expression can be written for Mg. — the twisting moment M r@ , and 


shear resultant vanish, this the formal solution for the elastic 
‘It may be noted in passing that when the parameter h/lg — 0 in ne. 23, the 
‘solution for the Poisson-Kirchhoff plate on elastic foundation is recovered, 
This form of solution has been previously given by Holl? and Woinowsky - 


_ VISCOELASTIC PLATE 0 ON N VISCOELASTIC } ‘FOUNDATION | 


linear viscoelastic behavior, the preceding analysis is extended. 

‘rect application of correspondence principle relating elastic and visco- 
elastic” boundary value problems leads immediately to the formal solution of 
Assuming | an operator-type stress-strain law for plate and foundation,9 
and following the correspondence principle, we replace the elastic coefficients 
— in all of the previous equations by the Laplace transforms of the 
emanate operators. At the same time the surface load, now a function of — 
time, as well as all other dependent variables are replaced by their Laplace 
transforms with respect to time. 
__ Denoting the Laplace transform of w(r, t) by W(r, 8), we recall 


under suitable restrictions on the function wel. t). D Denote the Laplace trans- 
forms of the plate and foundation time-dependent stiffness operators by ] D(s) 
and k(s), respectively, : and the transform of the plate parameter 8 by 8(s). 

“Then the correspondence principle leads from Eq. 15 to the following expres- = 
“Ala, s) [1 1+B(s) h (a 


Ks) +B(s) h? Ds) 
In a similar manner Eq. a and the Hankel inversion theorem om the Laplace 
of the foundation interface | 


“Ing, Archiv., Vol. 3, 1932, p. 250. 
yo “Rheology, Theory and Applications,” by em. Eirich, Vol. A. Academic Press, New 
York, 1956, p. 400, 


EM VISCOELASTIC PLATES 4 
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2 
a 
| {1 + B(s) a Ja J,(a@ r) 
«Thin Plates on Elastic Foundations,” by D. L. Holl, Proceedings, Fifth Internatl. 
Congress of Applied Mechanics, Cambridge, Mass., 1938, pp. 
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rie. a, ist the iterated Hankel- -Laplace transform of the surface 
pressure on the plate. | Similar expressions can be written for the shear re- 

- sultant and bending moments, replacing Eqs. 21, 22, and 23. In : actually per- 

- forming the inverse transformations to recover the time and space dependent 
variables, it is expedient to compute the inverse Laplace transform ‘before 
performing the Hankel inversion. Inasmuch as the Laplace transforms of the 
viscoelastic operators D(s), k(s), 8(s) are expressed as polynomials, the 
= inversion involves, in general, inversion of the ratio of rational func- 
tions. On the other = hand, performing the Hankel inversion first, if at all 


elasticity | by transform methods has been noted by Naghdi and W. C. Orth- 
- wein.10 In the following section an example will clarify the details of the 7 


‘plex plane. ‘This: feature ‘of the solution of axisymmetric problems in eco 


ve A CIRCULAR AREA 
As an example of the method of solution, the problem of an infinite plate 7 
Pet sr loaded over a circular area will be treated. Maxwell- -type visco-_ 
_ elastic properties will be assumed for plate and foundation. Furthermore, in 
the interest the plate and foundation taken ai as incom- 


pe no 


q H(t), 
| 


in which H(t) is the unit iterated Hankel- -1-Laplace t a 
pon 
of q(x, t) for this case q 


formed. 

a “Response of Shallow Viscoelastic Spherical Shells to Time-Dependent | — 
- metric Loads,” by P. M. Naghdi and W, C. Orthwein, Quarterly of Applied Mathematics, 
Vol. 18, Providence, R. July, 1960, PP. 10 107- 121, 
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re 
The terms Tp, 7¢ are relaxation times for the plate and foundation, respec-_ 
_ tively. The surface 
| 
ana into 20 and co and recalling tne definitions O 7 a 
Bo 
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manner, expressions” for the transverse ‘ear resultant 


moments can be obtained from Eqs. 22 and 23. 


has been carried out « on a , computer using FORTRAN programming fora num- 
ber of cases of interest. 


NUMERICAL EXAMPLE 


— 


‘The maximum deflection and maximum pressure 
was carried out for the preceding example by numerical integration, using 
Simpson’ 8 Rule, of Eqs. 30 and 33, with r=0. The geometric parameters i in 


the problem were to be Computations were for 


= , 5, , 10, for 1 values of = = 0, 0. 10, 1. 0, 10. 0. The results of these 


are shown in Figs. 3 4. first to Fig. 3, it is in- 
teresting to note that for the case 7,/7, = 0, corresponding to aviscoelastic 

_ plate on elastic foundation, the maximum deflection is quickly attained. 
the ratio 1.0, the response becomes largely linearly 


‘cous. In Fig. 4 it may be nc noted that for | the case P « = 1, the | pressure ¢ corre- 7 


— 
— 
| 
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FIG. 3. —CENTER DEFLECTION vs. ‘TIME FOR VISCOELASTIC 
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FI IG. 4, —MAXIMUM INTERFACE PRESSURE VS. TIME FOR VISCOELASTIC 7 
PLATE ON VISCOELASTIC FOUNDATION 
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"sponding to: an n elastic plate on elastic foundation is attained and maintained as 


plate foundation relax. Returning to the case ¢ of plate 
on elastic foundation = 0), a a limi 


the load were directly to the foundation. 
Another limiting case appears if Tp/T, = corresponding to an elastic 
_ plate on a viscoelastic foundation. The related problem of a thin plate (8 = 0) 


on a Winkler foundation (k = p/w) has been studied by B. C. Hoskin and E. H. — 


A of onde of bending for plates on a 
_ deformable foundation has been presented. The plate theory incorporates the - 
= of transverse shear deformation, and the possibility of viscoelastic ma- - 
properties in both plate and foundation is included. Using integral 
transform methods the deflection, interface pressure, and plate bending mo- 

at ments are expressible in terms of integrals containing geometric, viscoelas- 
tic, and dimensionless time parameters. This form of solution, in contrast to 
direct solution of the differential equations governing the system, is particu- a 
larly suited to numerical computation. For selected values of geometric vari- 
ables, a parameter study for viscoelastic and time variables is readily car- 
ried out. It should be pointed out that the integrals arising in computation of 
the bending moments in the plate are slowly convergent and asymptotic ex- 

_ Pansions of the Bessel functions occurring in the integrands is required. 
_ Finally it may be noted that a Maxwell-type viscoelastic model is unsuited — 


SUMMARY 4 


for describing the long-time behavior of a Physical = and is used here 
only for convenience of illustration. 
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APPENDIX. —NOTATION 


‘The following . symbols, adopted for. use in the pap paper, conform essentially 


with “American Standard Letter Symbols for Structural Analysis” (ASA Z110. q 
8- 1949) prepared by a committee of the American Standards Association with — 


‘Society representation, and approval by the Association in 1949: ai a 


modulus of elasticity 


‘a 


ce 


M,; 

= ressure at late- -foundation interfa 
radial t transverse shear stress resultant 


= transverse surface pressure on the plate 


= vertical nome of plate middle surface hay 


= Poisson’ s eratio 
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By Sciammarella ‘and A, J. Durelli 


—_ _ The objective of this paper i is to show the application of moire fringes in 
the determination of strains in two-dimensional models. The moire pattern is 
interpreted as a function of displacements only from which strains in the La- 
grangian and Eulerian descriptions are determined. This presentation is not 7 

- limited to small deformations, and sets of formulas appropriate for every case 7 


In recent years, moire fringes have — increasing application i in the field 
, stress analysis. Two are the main fields of application, deflections of plates” a 
and strains in two dimensional problems. In this paper attention is focused © 


from moire dataonthe latter, — 
_ In order to obtain strains, two different approaches have been followed, one 
that can be called “geometrical” and the other which consists in relating the 
_ fringes to the displacement field, The purpose of this paper is to present this _ 
“seometrical” interpretation seems to have been originated in a = 
- Published in Dutch in 1945 by D. Tollenard,3 and was applied for the first time 
Note,—Discussion open until July 1, 1961, To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. This paper is part 
7 of the copyrighted Journal of the Engineering Mechanics Division, Proceedings ¢ “ the 
American Society of Civil Engineers, Vol. 87, No. EM1, February, 1961, 
Assoc, Research Engr., Illinois Inst. of Technology, Chicago, Tl. 
2 Prof, Civ. Engrg., [linois Inst, of Technology, Chicago, Ml, 
Moire—Interferentieverschijnselen bij rasterdruk,” by D, Tollenaar, Amst 


> 


_ Instituut voor Grafische Techniek, 1945, 


— 
MOIR GES AS A MEANS ANALYSING STRAINS 
2 

— 
= 
7 
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ubject ‘of the strain determination in 1952 by ‘J. K. Kaczer and Kroupa* 
of the Physics Institute of the Charles University in Prague 
_ The basic idea behind the “geometrical approach” is the following. Moire > 
_ fringes are formed by two interfering line screens, One is printed inthe model — 
and is subjected to deformations produced by the applied loads, ' The second 
screen applied on top of the first is used as reference or master. . The moire 
pattern formation can be studied as the result of the intersections of the two 
s =_ mentioned systems of lines. Knowing the distance between the master 
_ grid lines and measuring the distance between fringes, _ it is possible by geo- 
metric a analysis of the intersections of the two _ systems of lines to. compute the 
distance between the model grid lines at a point, and the corresponding change 
in direction. With these two data, normal and shear strains can be computed. 
The “geometrical” approach | gives values of the strains that are the average 
- values in aregion limited by two fringes, because the applied formulas are val- 
Similar ap approaches are presented in other papers.‘ 5,6,7 Bromley8 obtains 
~ results similar tothose previously referred to by using! tensor notation. A dif- 
"ferent point of view in the analysis of the moire patterns was presented in 1948 
yy by R. Weller and W. Shepphard. 9 These. authors described the application of | 
x moire fringes asa means of measuring displacements. In 1954, M. Dantu,*” 
‘ Brounen the same lines, introduced the interpretation of the moire patterns 
=  interms of the components of the displacements, Dantu’s ; presentation i is lim- 
‘ ited to the field of small strains and deformations. A broader point of view is — 
presented herein, together wit! with the the corresponding set of formulas to b be used 


FUNDAMENTAL PROPERTY OF THE MOIRE FRINGES 
: q Two line screens are used in the strain anal sis, the “model grid?” ; and i the 
4 y “m g a 
- «i “master grid.” The distance b between the grid lin lines s is called pitch and 1 is rep- 

: ‘resented by yp. By it interference of the two screens a pattern of fringes is pro-— 
duced, The distance between fringes is called “fringe spacing” and has been 
assigned the symbol 5. Any line perpendicular to the master grid lines will 

: be. called principal section and a line parallel to the master grid lines, a sec- = 

ondary section. Model and master grids are assumed to be in the same plane, 


Let us consider the interference of two line screens parallel to each other ro 
be of different pitch, It can be assumed that one is obtained from the other by 7 


uniform contractionor elongation, This mechanism of fringe formation is shown. - 
. “The Determination of Strains” by Mechanical Interference,” by J. Kaczer r and w. 
Czechoslovak Journal of Physics, Vol. 1,1952,p.80. 
_ 5 “Untersuchungen zur Theorie der Doppelraster als Mittel zur Meszanzeige,” by R. 
Lehman and A, Wiemer, Feingeratetechnik Heft, 1953, pp. 5-199, 
+6 *The Measurement of Plane Strains } by a a Photoscreen Method, * by J. D. C. Crips, 
7 of Moire Effect to Measure Plastic Strains, by A. ‘Vinekier and R. . Dechaene, 
pr 8 “Two-Dimensional Strain Measurement by Moire, ” by R. Bromley, a 
Soc., Vol, 69, Part3B,1956,. 
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a 0 “Recherches Diverses d’Extensometrie et de Determination des Contraintes,” by _ aaah 
Mz, Dantu, Conference faite au GAMAC, le 22 fevrier, 1 


>. When two: opaque coincide, ‘there isa maxi- 


the center line of the first light fringe indicates a displacement in the di- 
rection of the principal section equal: to P if the center of fringe ii 

for the n th fringe, the displacement is np. It can be concluded that the moire 
_ pattern gives the relative displacement inthe direction of the principal section | 


to the final or | deformed shape of the body. The pattern indicates that 


--CENTER OF A DARK FRINGE 


= 


q the mechanics of continua, the description of the displacements with the final 


= 
coordinates of the points as independent variables is called euneetan descrip- . 


The so-called nominal reduction-area strain 


in Lois the initial length of ne specimen and L, , the final is con- 


sistent Eulerian 1 description, this definition tothe pattern of 
1 


engineering or con- 


oire data the engineering strain 7 


Fig. 1. A dark fring will appear 
— over a transparent stri + 
mumof light intensity a 
a: : _ grid is the one with the largest pitch, a point P in the undeformed state has u- 
center oF rence 
1.—FORMATION OF MOIRE FRINGES IN A TENSILE SPECIMEN 
= 
bas, one with the Lagrangian description isthe so-calle 


3p 


~ 
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PROJECTION OF THE GRADIENT 


wre 


hh 
| G THE SAME 
= j= RELATIVE DISPLACEMENT IN THE DIRECTION NORMAL TO THE MASTER © a 
= 

Kon 


x ee on easily shown that 


As shown in Fig. 2, the Moire fringes are the loci of the points presenting y 

the same relative displacement in the direction normal tothe master grid. Call — 
zero the first fringe to the right in Fig. 2, one the second fringe, two t the third © 

- fringe, a and so on, Assuming that the points of fringe zero do not « experience | 

- any displacement, the fringe one is produced by the displacement pof the points — 

+ of the model grid in the direction of the principal section, fringe two is pro-— 
duced by the displacement 2 p, , and the nth fringe by the displacement n p. . Then 
uf ‘the moire fringes are the lociof points with a relative displacement in the di- 
_ rection of the principal section which is equal to an integer number times the 
pitch of the master grid. These displacements are given with respect to the 
deformedor final shape of the model. Each fringe is characterized by a para- 
meter. This parameter is arbitrary since we are considering relative dis-_ 
"placements, This parameter is called “order of the fringes” and is assigned © 


In the following, a cartesian system of axes x and y is used as a reference 


symbol u and the component of | ‘the ‘displacement parallel to ery y- ~axis is as- 

Signedthesymbolve 

The component of displacement of a point in a two-dimensional continuous 

7 _ medium parallel toa reference direction is given by a function of two vari- 
ables o1(x, ye Here i can take the values 1 or 2, 1 if the reference direction : 


» 4 efer . The prior mentioned 
| has the geometric ‘The function z = $4(x, y) 

a in cartesian coordinates represents a surface (Fig. 3). | This surface can be 
represented by the projection nof it its contour 1 lines on the xX- -y plane. These lines 
_ a are intersection of the surface with planes of equation z=np, The e resulting — 


which k is a 2 constant. This is also the of the moire fringes. 


partial derivatives of i(x,y) are given 


in 


By 
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"FIG, 4.—CONSTRUCTION OF THE INTERSECTION CURVE OF THE | SURFACE 
WITH THE PLANE x = x, 


— 
CURVE $2 (Xo, ¥) - - 
PUY 


_MOIRE FRINGES 


ed To compute the sauteed derivatives of ¢; (x,y), the curves of intersection of 


_ the surface with planes of equations x = Cy and y = Co can be drawn, in which 7 A 
C; and Co are constants. The procedure used to obtain a partial derivative i . 

illustrated in Fig. 4, _ The horizontal position of each point of intersection of 
‘the moire fringes with the line AB (trace of the intersecting plane) is first 

_ projected on to the base line CD, The orders are read fromthe moire fringes 
_and distances equal to n p are scaledupf fromthe base line, A linedrawnthrough 
the points, thus plotted, defines the eens section, The | slope of this ‘curve at 


correspond tothe Eulerian There aretwo 
tions to the Lagrangian description. Assume that the patterns corresponding — 
_ two orthogonal directions have been determined, — A point P, of initial co- 


ordinates xp and ie! moves, after deformation, to point P with ‘coordinates x} 7 
and If the components of the experienced by the 


a: using the preceding relationships, ae is possible to obtain the necessary 
data to replot the moire patterns referred to the initial configuration, 
; a _ Another possibility is the following. The derivative of u with respect ‘to Xo 
can be  aeeaenes as a function of the derivative of u with respect to Xy he 


e 


“ee oeeeeeeeee 


"Likewise, it is to prove 


| 
9) 
ne 
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a Using the preceding equations, it is possible to go from one type of description 


STRAIN ANALYSIS 
_ _In the case of large strains, the deformations at a point is characterized in 
the description by the following and components of the 


7 In Eqs. 15, the subscript zero is ; used to indicate ‘that the initial coordinates | 
of the points are the independent variables, 
the Eulerian description, using the subscript to denote that the final co- 
ordinates are the independent variables 


‘From: the moire patterns in two J orthogonal directions, it is possible to com- - 
_ pute the derivatives appearing in Eqs. 15 and 16 by using the procedures « out-— 

_ lined previously under the heading ' “Fundamental Property of the Moire Fringes” z 
dl The components of the strain tensor given by Eqs, 15 or 16 can be given a 
physical significance which is is important when the classical theory of elasticity 
linearizations are made. In the following reference is ; made only to the La- 

_ grangian description used in the theory of elasticity, 2 
an Applying the definition of strain given in Eq. 3 to the elements of line MjNq | 

(Fig. 5) which before deformation are equal to MN * dig = 


MoPo = dy, =1, 


— 
— 
— 
4 
ee 
— 
: lle 8,” by I. S, Sokolnikoff, John Wiley and Sons Inc., bead ee — 


OIRE FRINGES: 3 


pe Fig. 5, it is also possible to shaun the sine oot the change in angle o 


Eqs. 17 can be applied to compute the normal and shearing s' strains from the 
7 data obtained fi from the moire patterns corresponding to two orthogonal direc direc- c— 


og 


G. 5.—D! 


eX 


‘Eqs. 17 and 18 can be xpan the in series, 


— 


4 4 


_ From Fig. 5, it is possible to see that the preceding simplification | is geo- 
‘metrically equivalent to replacing the strain of the line element M1Ny_ by the 
- strain of its projection on the x- axis. From Fig. 5 it is also possible to con- 


_ clude that if there is no rotationof the line element, tha that is if 8V0/8Xo sail 0, Eq. 


By expanding in se series Eq. 18 and manona terms higher than the first or- 


a Eqs, 20 and 21 are the oanatiann used in the classical theory of elasticity. 
Sometimes 8u,/8Xo is small and can be neglected in the series ody roe 


“of Eqs. 17a and b, and the shear strains are also small, but the rigid body ro- 
tations are large. This is found for example, in the case of the bending of a 
_-yery slender bar. The classical theory of elasticity can be applied to the prior 

= mentioned cases, but Eqs. 20 and 21 must be corrected for the effect of rota- 
tion, For moderately large second term in the series expansion 


With the same degree of Bq. 18 can he simplified t to 


— the range of application of Eqs. 20 and 21 or 22 and 23, from the 
Zs - moire data, it is possible in each particular c case to conclude whether or not 
the non-linear terms ‘should be retained, 


4 


4 SIGN CONVENTION 


This section is devoted to the fundamental rules of sign determination. 

By convention, the derivatives in the direction of the principal section are 
positive when the vienna of the projected line element is tensile. Tofind the 
Sign of the derivatives from a moire pattern, it is necessary to observe ‘the | 


_ Assuming that the principal section is inthe direction of the x-axis, the dis- 
_— ‘placement ofa a point at a distance L from the origin of coordinates will be 


f 8u,/8X,_ does not change sign at an intermediate point, the : sign of the dis- 


laceme nt andt the e sign 0 of the derivative are the same. From the continuity of 


— 
664 
unity, and 8v,/8Xo is so § wer Can be neglected with 
| 
— 
— 
— 
— 
| 
| 3) 

f 
— 

= 


the diateesiomnn on and their derivatives, it follows that, to change sign, the de- 
rivative must go through zero, Since, at every point the gradient is perpen- 
dicular to the contour line through that point (Fig. 3) and taking into consid 
- eration ‘Eq. 8, the tangent of the moire fringes is given et 


According to Eq. 2 25, when 


an 
te, 


ee 


{ablished To k know the signs, _itis necessary to know the absolute 
displacement of a point of the model, From Eq, 25, the sign of the derivatives © 
in the direction of the secondary sections is known, ‘if the sign of = giao 
in the direction of the principal section is known, 
vet 


ments by using sections parallel to the coordinate axes was ; developed. An- 
other possibility is to use the gradient and to apply Eqs. Ta and b, This method 7 


Inthe preceding, the method of computing the derivatives « of the displace- = 


‘—_ in the direction of the gradient can be used, If less precisionis required, 
= direction of the gradient can be assumed 1 to | be that of the « common normal 


—_ the common normal, and pis the = pitch, 


=a to the. type e of description ‘used because the following derivations are in- 


| 
7 


th Olre ge pent parallel tothe x-axis. Ifthe prin; 
——— cipal section is in the direction of the y-axis, and 8vQ/@ yo = 0, se 
and the moire fringes have a tangent parallel to the y-axis. 
Duetothe continuity of the displacements, all points corresponding to changes 
-. of sign are located on lines which are boundaries between regions of different _ FY 
. 
4 
4 
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Ing the following, the gradient at a point will be indicated by the —_— or. 


or ¢ n2 according to the case, The derivatives in the direction of the coordinate 
will be indicated | by 


‘These symbols repecsent, in general, the derivatives of the the displacements 
in the direction of the principal sections and the sum of the derivatives in the 
direction of the secondary sections. They represent the a of the 
strain tensor only when the strains are small. ‘aad: i 
_ The construction of Mohr’s circle is shown in Fig. 6. The. X-axis is as- 
sumed to be in the direction of the algebraicly greatest derivative. With 
convention, the angle between Oni me _ €- axis: ‘is the same a as the angle | be- 
tween 9,1 and the x-axis, 


Going from the physical plane tc to ie €-y plane, the y-axis ‘must be e rotated — 
clockwise a/2, ‘The vector rotates rigidly attached to the y-axis. 


Inthe e- -y plane, Fig. 6b, a pole Q is taken and, from this pole, the vector 
en parallel to 9,1 is drawn, From Q, the vector QM which is perpendicular 
to p92 is also drawn. - The di direction of QM is obtained from $n2 by a clock-— 
if the points M and N are joined, the vector MN is obtained, whose projec- 


tions > MP and NP are equal to ee a 


(Bla 


The ro jection NP ise equal 


=¢ m1 sin #1 n2 COs #2 = (38), 
which is at MP/2, is the center r of Mohr’s circle. The radius 


of the circle is obtained Dy joining C with The origin of coordi- 


"EXAMPLE 
As an example of the of the moire the of a disk 
7 Be diametral loading was performed, The solution of this problem using 
geometric approach was presented elsewhere, 1/2 in, thick sheet of 


12 “Geometry of Moire Fringes in Strain Analysis,” by S Ss. Morse, A. J. Durelli, and : 
A. Sciammarella, Proceedings, ASCE, Vol. 86, No. EM 4, 
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é 

FIG. 8 ION 


MOIRE FRINGES” 


a urethane rubber (Hysol 8705) was cast. A 300-lines per in, grid was photo- 
printed on the rubber using a photoengraver’s albumenoid base coating. After 
printing, a 4-in, disk and two tensile specimens were machined fromthe | 


ae: Two tensiletests were performed to find the elastic constants of the mater- 


ial. Fig. 7 shows the corresponding patterns at one stage of loading. In Fig. 
rsal displacement, 


=GRID LINES 


IG. 9.—DISK UNDER DIAMETRICAL COMPRESSION = ul 


DISPLACEMENTS 


in Fig. 7(b) fringes are the lociof points with equal longitudinal displacement. 


The patterns were obtained by the double exposure method, the unloaded and 
loaded conditions being successively exposed on the same film, valves 
_ The disk was loaded between flat plates, with the grid lines in the direction 


= 
OY 
— 7 
— 
4 of the load, Fig. 8, and with the lines perpendicular to the loading direction, — ‘ 
FFiig. 8 Moire fringes give horizontal displacements u, and in Fig. 
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ON. A TENSILE 
= 


EXPERIMENTAL 


EULERIAN ). 

LAGRANGIAN 
THEORETICAL 


nore’ AREA UNDER CURVI 


‘ 


RINCIPAL STRESSES ON THE “HORIZONTAL AXIS OF A_ 
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‘The loci of points: at which the derivatives in the direction of the principal 


and the secondary sections are equal to zero, are > indicated in Figs. 8 and 9 
ow On the horizontal diameter (Fig. 8), the displacements are positive and, there- 
- all of the central region enclosed by the line 8u/8x = 0 has positive de- 4 
rivatives. In Fig. 9 the displacements in the vertical diameter are negative 
and, therefore, the derivatives inthe direction of the principal section are neg- — 7 


ie sign changes. The lines of zero. derivative intersect the principal and sec- 7 
ondary sections at points called singular points. T\ Two types of singular points 
are indicated in Fig. 8. 
‘The data obtained from the tensile tests are plotted in Fig. 10, ii ‘this fig- ig 
‘ure, he results are expressed as functions of the Eulerian strainand the nom-— - 
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FIG. 12,—DIFFERENCE THE AL STRESSES IN THE HORIZONTAL AXIS 


inal stress. The strain is the so- called true 
stress-strainrelationship, 
In Fig. 11, the principal stresses alongthe horizontal diameter are shown in 
-dimensionlessform. The ratiox/ris usedas abscissa, in which r is the radius 
_ of the disk, The ratio of the stresses to the average stress Oayg = P/Dt is used 
as ordinates, The diameter of the disk and its thickness is denoted by D. _ ‘The — 
"experimental results are given in the Eulerian and | ‘Lagrangian descriptions, 
The Eulerian description was obtained directly from the moire patterns, i 
the elastic constants of Fig. 10 were employed to compute the stresses from 
the strains. The value of 0, ve was computed according to the initial dimen- © 
‘sionsofthe disk, 
ait The strains in the Lagrangian n description were obtained from the strains | 
in the Eulerian description by using Eq. 13, and the coordinates of the | points” 
were computed by applying Eq. 10. In order to obtain the stresses from the 
strains, values of the tangent elastic modulus of the stress-strain relation-— 


ship, nominal stress versus engineering strain, were computed, In Fig. 11, 
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_the abs-issae are referred to the undeformed sage of the disk, Also shown 
in the same figure are the results of the theory of elasticity. In Fig. 12, the 
difference of the principal stresses on the horizontal diameter, obtained by 
photoelasticity, is compared with the moire results. 

The maximum difference between the theoretical and the moire values of 0) 

on the horizontal axis is of the order of 4% of the maximum stress inthe disk, : 
a The same figure applies for 0}. The difference of 01 - 09, computed by the — 
4 moire method and photoelasticity is 1.5% of the maximum RT ns 7 


For a point with coordinates ink 


: 


a rectangular three- -gage, , 45° strain rosette analysis was also seutorened. mer 
this purpose, the moire pattern with a model grid at 45° with the line of load- 
a 
ing was determined (Fig. 13). The engineering stra strains at t the —— 
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ue 


“From Mohr’ s circle 
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The e rigid — ee of the point can be computed from the n moire data ae 


= 0.0262 sin 5° = 0,00228 

0509 sin 28° = 0.0239 
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ae 13,—MOIRE PATTERN OF A DISK UNDER DIAMETRICAL LOAD OBTAINED WITH a 


_ WITH A MODEL GRID AT AN ANGLE OF WwW THE LOAD DIRECTION 
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TION IN THE POINT x= 0, 2r,y=0,5r 

‘ 


. The rigid body rotation is equal t¢ to 
> -[0,0239-0.00228] = 0. 0108 


the of the principal derivatives 


to ‘Eg. 22, for the rotation of the element of length is 


his s quantity can in view of the pr precision in n this example, 
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‘SYNOPSIS 
Design procedures for minimum weight, statically determinate beams are 
using a deflection criterion. Optimum flange width and web thickness. 
variations are described for I-beams designed with respect to bending and 


shear deflection. The role of material ee beams is discussed 
from the point of view of stiffness and strength. ow Taye 


This paper is concerned with the minimum weight - -stiffness stat- 
“teally determinate bending members with specified depth variations. In 1936, 
Roxbee Cox treateda similar problem for torsional stiffness. 2 He concluded 
that maximum stiffness for a given weight is attained when the stress level is 
constant. In April, 1958, B. Saelman demonstrated that this condition generally — 
does not result in maximum stiffness for the torsion problem.3 Several monthe 
= to the Saelman article, a design procedure for optimum stiffness- -weight 

bending members was developed4 for Rock Island Arsenal, which procedure © 


- Note.—Discussion open until July 1, 1961, To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the | 
American Society of Civil Engineers, Vol. 87, No. EM 2, February, 1961, — Aen Se 
me... Assoc, Research Engr., Armour Research Foundation of Illinois Inst. of Tech, 
Stiffness of Thin Shells,” by H. Robes Cox, Aircraft Engineering, September 
3A Note on the ‘Optimum Distribution of Material in a Beam for Stiffness,” "by B 
Saelman, Journal of Aeronautical Sciences, April, 1958, p. 268, 
4 “Lightweight Structures and Prestressed Launcher Components, Rock Island Ar- 
senal, Ordnance Project No. _TU2- I Report,” — January, 
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: ; also demonstrated that a constant stress level is not generally a a sufficient con- 


os dition for maximum stiffness. This paper represents an extension of the latter 
a) Basically, the problem treated herein is that of adjusting one or two dimen- - 
_ sions of a beam cross section in such a way that the deflection at a given sta- c 
— 7” assumes a specified value and the total beam weight is minimized. wal 
_ The body of this paper consists of three sections. In the first of these, beam — 
be sections are treated for cases where only the bending deflection is considered | 
and only one dimension of the beam section is varied (one-parameter beam). 
For this problem, three distinct situations arise depending on the nature of the 
loading and the specified deflection. In one instance, no solution exists; in an-- 
other, a zero-weight beam can be approached which satisfies the specified de- 
flection; and in the last, the optimum beam dimensions are obtained by solving © 
the simplest problem of variational calculus. In the latter case, 


ed among their parameters— deflection, depth, span 
and material properties. For beams with large span- to-depth rat ratios, these a 


™  lationships closely approximate the ‘corresponding, but more cumbersome, re- 
. lationships obtained in the third section of this paper, which deals with beam | 
sections designed | for bending and shear deflection in which two dimensions of — 
the beam section are varied (two- wo-parameter beam). For certain n types of 5 sin- 
gle-parameter beam sections, the results developed in the first section can be 
modified slightly to take account of | both bending and shear deflection. | _ Example 
3 in the third section utilizes | such a a . modification fora ‘rectangular cross sec- 
In the second section of the body of the paper, constant- -strength beams are 
compared to the optimum deflection designs developed in the first section. 
i Conditions on the loading are established for which the two o types of woes are 


a ‘The third section considers three types of two- -parameter beam sections 
designed for bending and shear stiffness: (1) an I- beam withan optimum flange — 
width and an optimum 7 web thickness variation, (2) an I-beam with an optimum 
‘flange width variation and an optimum constant thickness web, and (3) an opti- - 
- mum prismatic I-beam. The discussions of the first two cases follow closely 
that given forthe one-parameter beam sections; however, the details are some- 


ie what more involved. In the third « case, , ordinary calculus is used to establish — 
_ the optimum flange width and web thickness for prismatic beams. The three _ 
_ types of deflection designs are then compared for a particular case. o 
_ For eachcase studied in the various sections, the method of design and per- 
tient formulas are indicated. A number of examples have been worked out to — 
demonstrate the design techniques: and to some of the 


‘Notation. —The letter symbols for use in this paper are 
they first appear, in the illustrations or in the text, and « are arranged alpha- 


betically, for convenience of reference, in the Appendix. 


"BENDING D DEFLECTIONS ¢ OF SYMMETRICAL BEAMS 


(SINGLE-PARAMETER BEAM SECTIONS 


which both minimises its weight and keeps the at anys station 
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‘BEAM DESIGN 


subjected to an arbitrary static loading in the plane of symmetry. It will be 
assumed that the bending deflection at a station x = £is given by the virtual- 


at a specified value. Consider a statically determinate, 


| = 
| 


> 


represents the modulus of elasticity, and isa variable ‘principal 
of inertia inthe plane of bending. 5 The integration extends over the entire span 

“ s. It is further assumed that the external moment M is is independent: of the the beam = 

- Three situations can be distinguished with regard to the ‘product M m: ( 1) 

] the sign of M m is not the same as the sign of the specified deflection A(t) ) 
anywhere in the span; (2) Mm assumes both positive and negative values in 
the values in the span; or (3) the sign of M m does not differ from the sign of 
ACE) anywhere in the span. We} shall consider each of these cases separately. . 

No physical distribution of material exists that will provide the specified 
deflection A(£). The satisfaction of of Eq. 1 can only | be achieved with a nonpos- 

a: ‘For this case, any specified finite deflection can be obtained witha alli 


of arbitrarily s small weight. The integral in ~ 1 can be: written as 


“positive 
which the first and second integrals extend portions of the span 
where the product M m is positive and negative respectively. An I' (x) and an 
-I''(x) are sought which will produce a ‘Specified 1 value of A(é). If I' (x) and © 
re ) retaken as constant in their respective spans, it is readily seenthat as _ 
4 I" is taken smaller andsmaller the magnitude of I' required to keep the right — 
‘side of Eq. at a specified positive value also becomes smaller and smaller. 
‘By reversing the roles of ¢ and I'' , the same argument holds w when a specified 7 


value is sought. any value of can be obtained with mo- 


‘further in subsequent sections. ot 
m/A(E)>O0 


stations other than x = increase 4 will be 


‘out loss in n generality, a as ‘non- negative 2 quantities. nt this case, it will be a 
sumed that the area and the moment of inertia of a beam section ba on one 


“York, ard Ed. p. 143, 


ie 
f 
= 
- ertia imply vanishing areas and, hence, an arbitrarily light bending member, . 
: 


in which A is the beam a area and p p denotes the weight density. This formulation 
_ q the simplest isoperimetric problem of variational calculus. 6 The stationary © 
character of the integral in 3is to the a 


x, t(x)] = 


are specified piec cewise continuous functions. ‘By substituting ‘Eqs. 6 and 7 into 


"i 5, one obtains an algebraic expression that can be solved explicitly for t(x) in 
terms of the parameter 


To cause this function to satisfy determined from the — ‘ob- 
tained Eqs. 6 and 


“Methods of Mathematical Physics, by and D. Hilbert, Inter- 
Publishers, , Inc., 1953, p. 216, 


. x parameter t(x), which must be determined. The beam weight, W, to be ua 
pa[t(x)) ax ...............0 
= 
— 


11) 
- Notice that t0(x) becomes a complex function when the product M m assumes 
both positive andnegative values in the span; hence, an extremal does not exist 
for this case. When M m2 0, the extremal does exist; however, it must be : 
shown that it minimizes the beam weight since the stationary character of the 
weight integral is not a sufficient condition for the occurrence of an extremum. 
From the definitions of terms in Eq. 10, it can be seen ~ 
bt(x)20- 


will be that w(t) 2 for that renders 


I= 
and satisfies the subsidiary condition A(&) = positive constant. 


Fe = W(t) - = 


Substitute for t(x) and clear terms: 

4! [a+bt(x)] dx - 


co 


ubstitute from 4 and 6 for. a an 


C(a+bt) dx - 


(a + bt) dx| 
° 


‘Since the quantities [M +bt)]1 1/2; and [ (e/b) (a + 1/2. are real, pos- 
itive, quantity in braces is non- negative 


ir 
4 
>, 
— 
— 
— 
’ 
~ 


a relationship between the beam weight w and the deflection A(E) ¢ can be 7 
obtained Eqs 11 into 3: 


vides the absolute 1 minimum weight. 


q 


sential can be casi in terms of weight by ee Eq. 18: 


oat 


| Because of the reciprocity of the isoperimetric problem, the the function tx) 
- given by Eq. 19 will minimize the deflection integral when the weight is spec- 
ified. 8 That is, the problem of minimizing the deflection with a given mass is a 
the same as the problem of minimizing the weight for a specified deflection. — 
Another easily handled form of the en of inertia and area is given by 


= 
are specified piecewise ‘continuous functions, n > 0, and k > 0. For these 
forms, the extremal function an¢ the weight- -deflection relationship are given 
by Eqs. 23 and 24, , respectively. Their development ‘par allels the previous case : 
precisely, but when demonstrating that the t°(x) given by Eq. 23 provides ab-- 
‘Solute 1 minimum weight, one should use Holder’ ~ Inequality? for integrals ra rath- 
4 


of Potential by 0. D. Kellogg, Dover ‘Publications, 1953, PP. 
we 


8 “Methods of Mathematical Physics,” _by . Courant and D, Hilbert, Vol. Inter- 
science Publishers, Inc., 1953, p. 258. 


9 “Theory of Functions of a Variable,’ by 

tions, Inc., 1957, p. 643. 
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when expressions more complicated than here are used 
for the area and moment of inertia, Eq. 5 cannot generally be solved explicitly 
for the extremal. Furthermore, a direct numerical solution of this equation ‘is 7 
not possible since it will contain the unknown constant were known, t 
could easily be determined numerically for ¢ every value of 2 . ‘This suggests — | 

d the following trial-and-error procedure: Assume a value of Ain Eq. 5 andcon- | 
struct numerically a function on t(x). ‘With this function, evaluate numerically 
the integral in the subsidiary condition, Eq. 4. If the value of the integral is ; 1 
equal to the specified constant, the function t(x) i is an extremal function; if not, 
as repeat the procedure with another assumed A. Since Eq. 5 usually has several _ 
‘roots for each value of x and several functions t(x) may be found for each A. 
_ Two examples will be treated to demonstrate the design procedure based co : 
Example 1: Ideal I-Section.—Consider the design of an ideal I-section, that P 
7 ‘is, two extremely thin parallel plates separated by a distance da. The area of ; 
each plate or flange at any x is therefore, 


~ AL Ally A 


sp 
deflection, and | material properties: of an in optimumly c designed ee I-section; : 


‘constant 


It can be seen trom 28 that the weight is to the de- 


optimum beams designed v with and titanium would 
all be equivalent since the specific stiffness is almost the same (E/p = 108 in 1) 
for most structural materials. This fact provides one of the motivations for — 


— 

eM 
he weight given by Eq. 2 
Note that the weight g 
° 
| 
| 
4 

— The optimum 4,ande=1; thus, 
11. where a = c = 0, b= d2 (x) /4, an 

O(x) = tO(x 2 a2] - 


considering ceramic demonstrate specific ‘stiff- 
nesses many times greater t than ordinary structural materials (boron carbide: 

E/p = = 7,22 x 108 in.). The optimum > flange area given by Eq. 27 is inversely — 
__ proportional to the modulus of elasticity; consequently, all other things being 
- equal, the lighter structural materials lead to bulkier flanges, which are ad- 
-vantageous from the point of view of stability. 


Example 2: Conventional I-Section or Box- Section.—Consider the m minimum 


weight-deflection design of a simply supported I-beam subjected to a downward-— 
— load system. The cross-sectional geometry : at any s Station is described 
in Fig. 1. The flange thickness r and the web thickness. ‘sare specified positive 
functions of x, the oe variation d(x) is specified such that d(x) > 2r, »and 
the flange width t(x) | is treated as the open parameter. The | moment of i inertia a, 
the area of section are given by Eqs. 6 and 7, 
. (29a) 
id 


4 


(29e) 


‘Fora downward- -acting load system on asimply supported beam, the external 

- moment diagram is continuous, convex with respect to the moment diagram 
baseline, and zero at the ends. Consequently, the product M m can always be 
taken as non- negative when a downward deflection is specified. The optimum 
flange width variation t°(x), then, is given by Eq. 11. 
As the ends of the b beam are approached, the product M m decreases con- 
tinuously to zero. . For this reason, Eq. 11 indicates that t0(x) is negative in 
some finite region at the beam ends. The stations where the flange width goes 7 
from positive to negative (cross — points) can easily be located by finding © 


The solution to a real aie demands, of course, that t(x) 20: conse- 

quently, Eq. 11 represents a ‘mathematical rather than a real solution tothe 

problem. If one formulates the original problem in a slightly different fashion, 

an optimum flange width va variation t*(x) may be 7 — to the additional 

For the beam section consideration, the of inertia and 
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‘Fig. 2 illu strates — a beam. The objective is to find an extremum t*(x) for 
_&<x<7n such that — and 7 locate its cross-over points. The expressions for 
and deflection of this" beam are by 


The minimization of to the condition that is equiva 
lent to minimizing X with Y = constant. This latter problem was solved at the 

_ beginning of this subsection; hence, the results { given by Eqs. 11, 18, ,and 19° 
- apply here if W is replaced by X, A(é) is replaced by Y, and the integration — 
_ over the span is contracted to the interval £ <x sn we extremal function 


for E< x <7. At the | over points, flange width is z zero; ‘consequently, 


and 7 may be determined from the conditions: 


)=0. 


(38b) 


(350) 


= 214.34 in.4 


— 
f 
4 
— 
— 
= 
> 
(34a) 
wae 
Numerical Example.—For a beam with the geometry and loading shown in 
3, the following quantities may be computed: 
— 


b, c, ande are constants. 


For these quantities and Eq. 33, conditions 34 Senn nal 
9 
(L- 


When this is ‘substituted into Eq. 30, that equation reduces to 


‘When this tte is for F P = 300 i in., 108 a 
214.34 in.4, and nd A( 150) = 


From Eq. 38, then 


is by Eq. llare 


_ EaA(t) _ 


Min, 


The attention given the stiffness and of constant-strength beams 
_ in the literature and the insight that such beams 3 provie~ about the properties of a 
m minimum - -deflection beams, furnish the motivauon for the the ensuing discussion. 


g 
(36) 
— 
The ratio of Eqs. 36 and 37 yields 
— 
q 
n= 
btained from the equation t°(x) = 0, where 
— 18.41 in, 
a ENGTH BEAMS 
CONSTANT 


1961 


(x) 


7 the neutral axis to the outermost fibers (q = d/2 > on By eames Eq. 40. 
} into the deflection equation, Eq. 1, the stress W% can be related to the deflection 3 


can be ‘to. the 


4 ) 


x, t(x)] = x) 

»~ resulting t(x) may then be used | to establisht the area AL: x, t(x)] and, hence, 


the weight. _ For the ideal I- section, the moment of inertia given by Eq. 26 is 
set equal to 1°(x). The open parameter then becomes 


| Mam, 


To clarify comparison of the weights of the o optimum “deflection beam and the 
constant strength beam, the weight of the optimum deflection n beam, Wa, given 
| 
ae “Strength : of Materials,” by S. Timoshenko, , Vol. I, Dz Van Nostrand Co., Inc., 


Eq. 28 28 may be rewritten as 
1955, p, 210. 


For a constant-strength I © cross section 
a loaded in a plane of symmetry, the moment-of-inertia distribution is given by Ps 
my «the expression, | 
4 a 
using Eq. 41, the moment of 
a 
= For beam sections that are defined by one of —— 
if 
g 
ht 
and the weight of the constant strength beam Eq. 25, becomes > 
ae 
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. 4.—COMPARISONS AMONG PRISMATIC, CONSTANT STRENGTH, AND 
DEFLECTION FOR A SPECIFIED TIP DE- 


In Table of prismatic, constant strength, and minimum de- 


EML 
| 
a are proportional. Thus, the equality holds only when M/m = constant. 
if 
— 
Deallis are compared respect to 1oacines alice! - 
a tions. For the more common loadings on simply supported beams, it can be ae 


those of the n minimum - -deflection In ‘these cases, one would e ex- 
_ pect the properties of the minimum- -deflection beam to be somewhat similar to 
a _ those of the constant- -strength beam; that is, the stress level should be nearly — 
constant, and the deflection curve of a uniform- -depth member should be nearly 
circular with the maximum deflection occurring at the center. The cantilevers 
in Table | 1 show considerable disparity, however, between the weights of the 
constant "strength and minimum deflection beams. As is illustrated in Fig. 4, 
‘the differences become more pronounced with order of loading. 
-BENDING AND SHEAR DEFLECTIONS OF SYMMETRICAL BEAMS 
ct (TWO-PARAMETER BEAM SECTIONS) 
-_ In the previous ‘material, only the bending ; deflection was considered. ‘Here, 
_ account is taken of both bending and shear deflections | and three cases are 
_ treated: (1) an I-beam with an optimum flange width and an optimum ° web_ 
thickness variation, (2) an I-beam with an optimum flange width variation and 
an optimum constant thickness web, and (3) an optimum prismatic beam. 
‘The shear deflections of a beam _ be approximated by the virtual- -work 


in which Vis the extrenal ishear, v v denotes the shear due i a onatting unit load 
at station : x = =§,Gi is the modulus of rigidity, and a dimen- 


and s, t, r, and d are the dimensions ns shown in in Fig. ay 


Case 1: Optimum Flange Width and Web Thickness Distributions. _The mo- 
ment of inertia and the area the I- beam section ‘shown in ‘Fig. lare, respec- 


tively, 


A(x), = t(x) +h 


_ “Deflections of Beams with Special Reference to Shear Deformation,” by J, a. 5 
in and G, W. Nat. Adv. Comm, hare., 180, 1924, 


4 
sections, a = 6/5; for I-sections or box-sections, a is given 
— 
i Ve. (d-r)(d-2r)>0 ........... (50) 


‘an the use of Eqs. 52 and 53 and the form factor neni by Eq. 4, the beam 
weight anc and the beam become i 


“be noted that the unit load must be applied inthe same direction as the specified 
/ deflection. if the dimensions t and s are required to be non-negative, Eq. 55 7 
- reveals five distinct possibilities with regard to the products M m and V ve 


Eq. 55 may be written as 
— 


Regardless of the value of M m/E, the ratio t/s may be set equal toa constant _ 


of sufficient | magnitude to cause. the numerator of the integrand of Eq. 56 to 
change signs in S. With the dimension s still at our disposal, we have a situ- 
ation similar to tothat discussed in Case 2 under the heading “Bending Deflections — > 
of Symmetrical Beams.” ’ Consequently, any deflection can be obtained with a 
beam of vanishing lateral dimensions. = 
Case 1-b.—Vv<0 throughour S; M m >-Vvt E/G somewhere in 
Clearly, i in any portion of the : span, , the ratio t/s | can be set sufficiently large. 
to ensure that the numerator of the integrand of Eq. 56 is negative. If Mm> - 
; . v f E/G anywhere inthe span, the numerator can be made positive by setting 
t/s a: 0. Since the numerator can be m made p positive in some portion of the s span q 
bi aagennee in other portions, the circumstances are equivalent tothat of Case 
a. Thus, a member of vanishing weight ma may be obtained | which will ane 
Case 1-c.—-V v<0 andMm<-VvfE/G throughout S. = 
a. the conditions defining this case are fulfilled, the numerator of the inte- 
grand of Eq. 56 will be nonpositive throughout the span regardless of the value 
_ assumed by the ratio t/s. For this condition, the integral in Eq. 56 can yield 
§ a positive value for the specified deflection only if negative values of t and s 4 : 
admitted. Since this violates the ‘original hypothesis, no solution exists 


‘for this case. 


 - 
— A( E) “EI dx bt+gs dx ce 55) = 
functions, t° and sO,are sought that will minimize the weight W while sat- 
= @6=6=-—ssiinstead of dealing with the two open parameters t and s, we shall consider _ .* 
_the equivalent, but more convenient, parameters s and t/s in which s > 0 and 
> 
4 
— 
ACE) = : 
— 
as 
a 


Case d.—V v> 0 throughout § S; Mi m Vvf E/G somew somewhere in S. 
Inthis case, the ratio o t/s can be set sufficiently great to ensure that the 
merator of the integrand of Eq. 56 is positive in any portion of the span. | _— 
thermore, if t/s is set equal to zero, , this | numerator will be negative some- 
where in the : span. Again the circumstances of Case 1- -a prevail and a beam - 
vanishing weight that can achieve any | Specified positive deflection can be ap- 
"Case l-e.-V v>0andMm2>-V vfE/G throughout S. 
“4 The conditions defining this case insure that the numerator of the » integrand» 
of Eq. 56 is non-negative. | In this case, , it is possible to find two non-n gative 
functions t(x) and s(x) such that Eq. 56 will be satisfied. We seek those func- 
tions which, in addition, will minimize the weight integral, Eq. 54. This prob- 7 
lem is an isoperimetric problem in two arguments. _ The two extremal func- 
tions t0(x) ands0(x) must be determined to satisfy the following reduced forms 7 
of the Euler equations for a suitable value of the sapremnnirad Re. _ This eed | 


— tt + 


M m G 
‘Vv 


(&) 


; 
q 
| 
7 and 58 andsolving the resulting 
T 55; hence, 
— dx . .(61) 
— 
T variations 
ACE) 'G(hb - ge) E b | 
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_ With the use use of Eqs. | 62 and 63, t the optimum beam weight becomes 


E/G =25 5 for most ‘structural materials, weight is inve 


to specific stiffness, E/p. 
3 _ Noting that the quantities (hb - ge) and (fb - jg) are positive when the | 
depth d is greater that 2 r, one can readily show that the extremals cannot be — 
_ realandnon-negative when any of the conditions given for Cases 1-a to 1-d are 
fulfilled. Examination of Eq. 62 indicates that the web thickness s° is always 
- non- negative. On the other hand, Eq. 63 indicates that the flange width t° is 


non- non-negative only pwhen 


im 2—G 


it is possible, as was demonstrated in Example 2, to define the beam geometry 
in those regions where t° would become negative. In Example 2, the flange — 
_ width was allowed to decrease to zero; here, a slightly greater weight savings 7 
can be realized by permitting the flange width to decrease no smaller than the 
web thickness (a rectangular section is more efficient than a degenerate I-sec- 
tion). _ Thus, inthose regions where t° would normally be less than s°, the sec- - 
tions will be taken as rectangles « of width s(x) and depth d(x). The design p Pro- 
_cedure will be illustrated in the following example. 


Example 3: Simply Supported . I- “Section or Box -Section. a simply 4 


s° in the regions adjacent to the similar tothat 
in ‘Fig. 2, regions of sucha beam are > identified: region 0 <x< 


_ tion of each region will be denoted by Wi and Aj respectively, where the sub 


peript identifies the The total beam weight is 


wew, 
7 


~ tion are distributed among these regions. “The weight and deflection cee 


| 
— G(hb - ge) G \ — 
rselv pro- 
7 — 
— 


(67a) 


t(x) = s(x 


t(x) = (x)].. .(67e) 


nSx<L. 


a to the conditions Ai = = Cj. . The extremals or 3 for the regions 1 and 3 are 
obtained from Eq. 23, when the following replacements are _made: t0(x) “by 


for 


Wi = 0<xx [note 
; 
— 
90x); and S by 


region 2, the extremal $5 obtained from Eqs. 62 and 63, 
respectively, where S is replaced by x s 
| 
Gb 


“inwhich 


il Ea. 70 through 75 ‘completely « describe the beam geometry wcities the quanti- 
‘ties Ay, Ag, Ag, —, and n are specified. The values of these quantities are de- 
termined from the expressions for the total deflection (Eq. 68), the continuity 

of the web thickness at stations gE and n, and the equality of the flange and web 

_ widths at stations and Ne Hence, 


#8 


da, 


3 te 
83 (x, 43, ")" dL b+g VL beg 
— 


BEAM. 


‘The following numerical example for the conditions — in ses Swi a illus- 


Example.—From symmetry, Ay= 43 and the use 
of Eq 78, 42). Substituting for for M M, v, and v in Eqs. 70 through 


When 


‘With this of bo is found from E ‘79 to be 
— 


— 
N=— VE log E(f+j) E(f+)) +1 (85) 
eee one substitutes Eqs. 86 and 87 into 81 and clears terms, £ becomes a 


EM 1 


, the beam defined Eqs. 83 


through 90 is is plotted in. in Fig. 6, 6, with E = = 30x psi, G = 12x 108 p 
of 2 of this beam can be computed from Eq. when 


(a Loading Section Geometry 


= 
FIG. 5. _—MINIMUM- WEIGHT DESIGN OF A CONSTANT- -DEPTH I-BEAM_ 
is replaced | b < x With the use of the replacements made in the com- 


of 3 Eq. 24 gives the weights of ‘regions 1 and 3. The pertinent 


for this beam are: 


283 lb per cu in., 


= %3.527 (10)-5 x2 + 
= si, 
= 0. 05459 in. 
2.833 (10)- + 5.557 (10)-3 - - 0.2374 in, 


= 178.5 543 
— — 


= 0. 0.002995 in. 


4940 in 


= 0.5000 in. 


= 0. 5000 
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2: Optimum Flange- Width -Constent 
“Web.—The difficulties inherent inthe manufacture of a variable-thickness web | 
“plate lead naturally to an investigation of constant thickness webs. Once again 
_ for the I-section of Fig. 1, the moment of inertia andareaare =~ 


| 
in atte: s is taken as ; constant and 'b, g, e, and h are defined in Eq. 29. With 
; _ these stipulations, the beam weight and deflection are given by Eqs. 54 and 55, : 
respectively. In this section, we shall consider only those cases where an ex- 
= tremal function t°(x) exists. Proceeding formally, the condition for minimizing 
the weight W when A(£) is 2 specified positive constant is given by Eq. 57 alone. a 
_ When one performs the operations in Eq. 57, solves for t, and eliminates by 


i fim [y me, VveK, 
Ss 


Ebe’ Ge L G ™ 


Itis clear from Eq. 93 that the extremal fun function exists if, 


)= f Gb 


s. Note that t°(x) is always negative if 


_ Physically this means that the web is so flexible in shear’ that the total deflec-_ 
can never be to A(é). Even 


— 
% — 

(93) 
tm a 
ndif 


it is possible for £(x) to some regions of a beam. in 

the previous cases discussed, these regions may be isolated and their flange 


onstrated for a cantilever I- -beam under a triangular loading. ‘The flange | ane 
is restricted to be or to the constant web thickness s. 


ucts Mm and Vv increase zero at the tipto a maximum value at the sup- 


port. Under these circumstances, Eq. 93 indicates that t9(x) is negative in | 
some interval at the ‘tip, ‘say O<x <é. In this interval, the seas section is 


4 


(a) Loading (b ) Cross Section 
FIG, 7.—CANTILEVER I- SECTION UNDER A TRIANGULAR LOADING 7 


. 


me = 


placed by Y and the interval S is ‘contracted to0<x< E. "The weight of the re- 


— 
— 
q 
| 
i = 100 kios 4 = 
— 
x 
ae _ taken as rectangular, the width as s, and the depth as d. The weight andde- — i 
of this beam are given by 
ry 

cd :. AS in Example 2, the minimization of W subject to the condition A(é) =con- § 
is equivalent to the minimization of X subject tothe condition Y = constant. 
ad wih 


Mm ,Vv 


0 


‘Consequently, 


‘from which £ can be determined. Fora given s, Eq. 103 ‘must, in 1 general, be 

_ The optimum s may be found by a straightforward, mou cumbersome, 
procedure. Because Eq. 103 canbe solvedexplicitly for s, this quantity is sub- 
stituted into Eq. 101 yielding an equation for W* as a function of the independent , 
variable &. The quantity | W* isthen minimized by differentiating it withrespect 
to g, setting the result equal to zero, and _— that equation for the —— i 


a yields the optimum 


7 pre ‘ented for the cantilever beam shown in Fig.7. It was s derived by guessing 
ava a ‘ue of £; substituting it into Eq. 103 and solving explicitly for s; computing 
zg es u. ‘ing the resulting s anc s and | the original; and repeating the procedure with 


Numerical Example.—Using the loading and geometry shown in Fig. 7,a tip 
- deflection of 1/2., and the same values of E, G, and pused in Example 3, Eq. 


sulting beam can be found from Eq. 98 when W° is rep 
by Y,, and S is replaced by the interval O< x < & 
— 8 * 
ae, Summarizing, the optimum beam is defined as - — 

t(x) + (102a) 

Sie 
— 
ay: 
— 
4 
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BEAM DESIGN 


“7 12(b + g) ACE 


[assou E V2. 72: 432 14815 5 - 2. 179 18 


12, E+ 6032093 
ti 


(0. 21130 1, 2432) 
V2.1130 + 12, 12.432 
where W* is expressed in pounds ; and s s and E | in inches . By using Eqs. 104 aa 
105, values of s and W* were computed for different values of . The resulting — 
relationships between w* and ~and between w* and s are shown in Fig. 8. | a _ 
is apparent from this figure that the weight W* increases very rapidly as s is 
decreased from its optimum value or extreme point. On the other hand, for. 
ee values of s greater than the optimum value, , the curves are very flat. This in- 
sensitivity is very useful in the sense that the web thickness can be increased 
for strength payee with little increase in the weight over en 


found from Fig. 8. ‘It is s characterized by = 227.205 Tb, Simin = 0,.1800in., =f 


| 


al 
103 is specialized 62 q 
— 
where s and are expre: - 
J, G ) ....(105a) =o 
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in = 22 8 in, When the values of Smin and Emin are 8 substituted into Eq. 


t(x tmin = * V 


t(x) = 


— 


- Case 3: Optimum Prismatic Beam. -The cross- sectional geometry of the 

prismatic beam considered in this subsection defined in ‘Fig. where r and 
dare specified constant dimensions and s andt are open constant parameters. 


_ The moment of inertia and are of this beam are a etm 


where x and tare in inches. The plan and elevation views of this beam | are” 


24 


(106b) 


The weight and expressions become 


. (108) 


| 
lal ’ 


Our objective is to find values s and t ‘that will minimize the weight WwW sub- 
ject to the subsidiary condition that it ACE) be a specified positive constant. . 
one applies Lagrange’s method of undetermined multipliers, the optimum | 
mensions s° and must 108 and the following 


(t 


and Integral Calculus,” by R, Courant, Vol, II, Interscience Publish- 


— 
- 
‘ 7 
— 
and 
te 
ne — 
im 
il 
— 
4 ae 


* Optimum Flange ‘on Thickness Web; Weight 22721 Ib 
ex 84790107 - 1.3689(10) *x; 9.1856 $x 100 
> 


es 


10, MINIMUM WEIGHT ¢ DESIGN D DEFLECTION - 
DEPTH - 15 IN., FLANGE THICKNESS - 1/2 IN., LOAD -100 kips, SPAN 


— 


106 4 ‘February, 1961 


After performing the operations in n Eqs. 111 a d 112 and eliminating Lagrange’ 
multiplier A with Eq. 108, the optimum semaine s° and t° and the associated 


- When these results are ‘ne to the beam treated in the numerical example 
_ Of Example 3,s° = 0.06157 in., t© = 3.20676 in., and WO = 229. 880 lb. Recalling 
that the beam of Example 3 weighed 179. 624 Ib, we find that the optimumly tap- 


Fic and optimum constant thickness web, and F optimum flange and opti- 
mum n web. Case ) represents a 40. 6% weight savings compared t to Case _ 


The e investigations here herein have dealt with beams ns of specified 
epths and | flange thicknesses. For the elementary deflection theory —— 


_ here, one would find that the optimum flange thickness is vanishingly small 
and that the optimum depth is infinitely large. _ Changes in the flange thickness | 
have very little effect on the beam weight as long as the (d/r) - -ratio is rea-— 
‘sonably large; however, the weight is strongly influenced by the depth. Eq. 28° 

indicates that the weight of a constant-depth I-beam is inversely proportional _ 
_ to the square of the > depth. It should be pointed out that single-parameter beam 

- sections exist that can be proportioned by using the techniques developed with- 

; 2 out specifying the depth. For example, the optimum radius variation to(x) of 


solid circular beam can be found from Eq. 23 whenw=7/4andz=m™ 


_ Frequently, the « optimum web thickness based on a deflection design is so 
small that web buckling is experienced. Since the specific stiffnesses of most _ 
a common metals are the same, the relationships given herein indicate that ost 

5 use of alower density metal will give a thicker, and consequently more stable, 
design without an increase in weight. But, the flange widthvery often 
_intolerably large for such a material substitution. . In these situations, it may 


be possible to use a | bimetallic beam, with : a high- -density flange ‘material ‘and 
low density web materi 


2 
“palé)| G VGI\Eb G/) ..... (113) 
beam. In Fig. 10, the weights and geometries of the cantilever beam considered 
mm 44.4% (y= 3). Case (b) provides a 36.4% weight savings compared toCase 
= (a), only 4.2% less than the optimum flange-optimum web combination. _ a — 
— 
— 
| 


‘In general, if an element of a a beam | is replaced by an | equal- depth inal 


of a different material, the bending stiffness will remain unchanged if the width 
of the new element is given 
tg E 0/Po Po 
where the subscript n refers to the new element and the subscript o > refers to to 
the original element. When the specific stiffnesses equal, 


which indicates that the weight is for a bimetallic member of equi-— 

se stiffness. The same argument holds for the shear stiffness since the 

(E/G) -ratios of most common metals are the same. hana Pee es 
_ Several investigators have studied the design of constant strength beams 
wherein the weight of the beam itself was considered part of the loading. Ih 
this paper, the weight of the beam is neglected; consequently, the results are 
" exact only when the live-load deflection is specified. However, the deflections - 
fue to the beam weight are usually very small compared to the live-load de- 
flections and can be minimized or a by cambering or prestressing - 
As indicated in Example 2, the results described for I-beams apply without — 

- modification to box beams, _ It is often advantageous, from the point of view of _ 


_ lateral buckling, torsional stiffness, and shear lag, to use such multiple- web 


writer would like assistance of Miss» Anne Hum- 
_ phreys, Technical Assistant, Armour Research Foundation, who carried out 


= - specified f function; in the I- -beam, the moment of inertia of the web; Eq. a 


ha = area under r the Vv diagram; “defined b; by ‘Eq. 110; 
tue the I- -beam Eq. 29; 


= = specified function; in the I-beam, the web ; area; sled 29; also subscript 


“Strength of Materials,” "by s. Timoshenko, Vol. 1 I, D. Van Nostrand Co., inc., "3rd 
1955, p, 217, 


E REAM D N <> 0 
(117) 
— 
| 
g 
| 
- 
_a 
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February, 


= area under the Mm diagram; by Eq. 1. 109; 

= total depth at any section of an I- ‘beam; also aw denoting deflec- 


tion design; — 


of 


pages specified function; 1 in the I- beam, twice the flange thickness; Eq. 29; 


defined by Eq. $1; 


= modulus of 3 


de ined by Eq. 


“used as asap power in Eq. 20; 


moment; 


| = defined by Eq. 85; 


t denoti lel t: 
script den enoting original element element; 


subscript denoting prismatic c beam; 


= distance from neutral axis to outermost fiber; 

= total beam 

od 


2 


flange w width of an I-be beam; ope open parameter; oy alla > 

external shear; 


— 
ah, 
constant, used aS a power in Eq. 21; also subscript denoting new ele- 
oe 0 = superscript denoting association with an extremal function; also sub- 7 a 
on or open 


- specified function; us used in Eq. 20; . 
= defined where used; see Eq. (31 or Eq. 98 98; 


coordinate » along a a beam; 


= superscript denoting association with a modified extremal function; 


= "define where used; 

‘= defined where used; > 


tector; defined by Eq. 49; 


= cross-over to a | 
specified either deflection or tot 


shear deflection; 
‘specified positive stress level; 


constant or Lagrange’ s multiplier 


weight t density (weight /unit volume) ; end 


degree of the bending moment in Fig. 4, 


a EM BEAM DESIGN lo 
-Y  =_defined where used; see Eq. 32 or Eg. 99; 
— 
i: al deflection; 
— 
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& STRESS ANALYSIS OF SHELLS 


merical examples are Pot 


pe _ This paper i is concerned with the problem of structural analysis of thin ie 
_low shells in the form of a translational surface. In recent years, such shapes 
- the hyperbolic and elliptic paraboloids have beenfrequently used inthe field 
of thin-shell construction, and much attention has been devoted to = study of 
the structural behavior of such shells. 
The problem of bending of this « class of shell surfaces can be studied with 
e the scope of Marguerre’s general theory for shallow shells.? In this paper ‘ 
_ the Marguerre equations of the small deflection theory are specialized to shells 
in the form of a parabolic translational surface, resulting in a system of two 7 
fourth order partial differential equations with constant coefficients. The Mar- 


“ee —Discussion open until July 1, 1961. To extend the closing date one month, a 
4 pow request must be filed with the Executive —— This paper is part 


Struct. Engr., Oslo, Norway, Presently (1960) research at Univ, 
a California, Berkeley, Calif. 


2 “Zur Theorie der gekradmmten Platte _grosser Form&nderung,” by K, 


EM1 @G 
> 
nee _ A method of solution of the problem of bending of thin shallow translational 5 
— shells is presented. A generalized Levy-type solution including hyperbolic and 

q 
q 
q 


‘tern systems of three simultaneous differential equations for the displace 


ment components u, v, and w w resulting from the two formulations will be given, — ; “4 a 
from which it will be apparent that the Flugge- ‘Conrad equations are most con-— 
venient i in he salaagay apne and further developments in this paper will be based on 


from which all the remaining unknowns of the problem of bending can be « ob- q 
tained by differentiation. It is shown that the stress- displacement function has" 
_ to satisfy an eighth order partial | differential equation, This equation is inte- 


grated using the Levy- -type solution, applicable to a class of boundary © value 


problems. _ By inclusion of a shape parameter, the solution is valid for — 
_ bolic and eliptic paraboloids and cylindrical surfaces. 


-_ acteristic coefficients for the statical quantities in the shell in the manner of 
~L. Holand, 4 by which the computational work in application is greatly reduced. ; 
‘The form of the solutionis such that it lends itself totabulation as in Riidiger- 5 
‘Urban, 5 where influence coefficients are presented, The solutions are well 
suited for comparative shell shape studies, Tables of the roots of the char- 


acteristic equation are included, | The method is illustrated with examples of © 
uniformly loaded shells with various values of the shell parameters.  —~ 


It should be noted that the theory used is a shallow shell theory, which is | 
on the that squares of the first order 


derivatives of the 

surface functions, (z, x)? » (z, y)* , and (z,x)° (z, y)» are negligible in compari- 

gon with unity. This limitation on the application of the shallow shell theory _ 
_ has been examined by E, Reissner, 


who states as a general rule that shallow _ 
~ shell theory will be more than sufficiently accurate as long as (z,x), 


< 
(ay) 
1/8, and ¢ often accurate enough | for practical purposes as long : as (z x) (Zy )s 


_ Notation. —The letter s symbols adopted for use in this pa paper are defined where 


~ they first appear, in the illustrations or in the text, and are re arranged a. alphabet - 

ically, for convenience of reference, in the Appendix, 

OF THE PROBLEM 

Geometry of the Shell: Differential Equations. Consider in the following a a 
class of thin shallow translational shells, rectangular in aad the middle sur- 
face of which is in cartesian 


on the sign of Eq. 4 represents a a hyperbolic paraboloid > 0), a parabolic 
cylinder (y= 0), or an elliptic paraboloid (y <0). 


3 On the Calculation of Shallow Shells,” by W. Fligge and Conrad, ‘Stanford 
Univ., Div, of Engrg. Mech., Tech, Report No. 101, 1956. eB a! 


n 
_ 4 “Design of Circular Cylindrical Shells” by I, Holand, Oslo Univ. Press, Oslo, 1957. 


coordinate planes leads to inconvenient expressions for the in-plane displace- 
ments W_ DA Conrad? arrived ore able exnressions by 
f 
a 
— 
| — 
— { 
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EM TRANSLATIONAL | SHE LLS 


small deflections and a a ‘distributed load y), ‘the Mar- 


in which w is the maton wanes in the z- iit and F is an Airy 


= 


"Substitution oft the surface Eq. 1 into = 2 cet 


introduction of a stress- displacement function (x,y) defined by 


system of Eqs. 4 
qs. reduces to the single equation ; 


- All the unknown stresses and displacements inthe shell are given a term oe 


the stress function (x,y) by the following relations: 


-D( Waxy = 


— 
—- 
— 
| 
a 
— 
: 
which | 
(6) ff 
— — 

— 
&g 


1961 

=D 
+ v4(0, 


where the notation of Reissner® is employed; in which v ‘refers to Silents - om 
ratio;  Mxx, Myy are flexural moments; Nxx, Nyy refer to stress resultants’ 
- tangential tothe shell surface; and Nxy is the shear stress resultant tangential : 

- to the shell surface. It should be nen that to the degree of approximation in- 
oe in the use of the shallow shell theory the difference between tangential — 
stress resultants and stress resultants in the xy-plane is negligible. Stress” 
resultants in the z-direction, however, must incorporate both the transverse. 7 
_ shear stress resultants and the components in the z-directionof the tangential 
stress resultants. In Eqs. 7 the transverse shear stress” resultants and the 
oe effective edge stress resultants are denoted a. Qy and Rx, Ry re 


_ Spectively, whereas the corresponding stress resultants in the z-direction are | 


denoted Vx, Vyand Vx, 
Ss displacement components u and v are determined from the following | 
system of displacement equations which has been derived | from the system of 
equations given | by Fliigge-Conrad; 3 is 
2 


+2V2, xx2, yy*Z, yy’ 249(1- v)z,. 2) w-(h2/12) viw w = = q, (1- 


relations between u, v and w for qz = 0 
é 


a 


q = Q, + Mx, (Ti) 
= Qy May (7) 
| 
— 
r _ which is equivalent to the system given in Eqs. 2. Introduction of the surface ae 
aa function given in Eq. 1 and elimination of u and v results ina signel equation _ Lage 
mination yields the 
7 
= 


TRANSLATIONAL SHELLS» 
k 2 (- (1+y) (9 2w ), y- (1tv) (9 


which by introduction of the stress- displacement function, Eqs. 5 yield 
ky 2 (v2 9) C40) ( 
op & 
ve kp | ((1+y) (v2), 6) (10b) 
~The displacement equations resulting from the Marguerre formulation in 
terms of projected displacements U, V and | WwW corresponding tc to th the > system of of 
appear inthe form 


+1 2) (1-v) 1/2) 


in which capital letters have been used to distinguish the projected apfene- 
ments from the normal and tangential components. 
7 | timination of U and V from the system of Eqs. 11 again yields the a. 
ential Eq. 6a. The resulting components U Vv, however, 


given by the relations Pate 


which may be verifies by geometric consideration, Inthe ‘solution of 
boundary value problems Eqs. 10 will be used, 


i 
ph, woe 
— 
Which results in the foll 
i 


_ Thus, the problem of bending of thin shallow translational shells with ‘i 


tributed ' surface load of intensity az is reduced to the solution of the partial © 
_ differential Eq. 6a subject to appropriate boundary ‘conditions on the edges. — 
Boundary Conditions... each edge of the shell there are four independent 
_ boundary conditions, two of which are stated in terms of the in-plane stresses 7. 
and displacements and two are stated interms of the flexural quantities of “a a 
In the following it is assumed that the two opposite edges x = + a/2 have ee 
ment free support, and that the edge stiffening members are rigid in the di- _ 
rection of their axes and have negligible bending real resistance | in planes: tangent 


which, interms of the stress- function yields 


In the following an an edge having the e given in Eqs. 14(a) will be 
denoteda simply supported edge, 
™ solution of the differential equation is sought in such aform that arbitrary 
boundary conditions can be satisfied along the remaining edges 
SOLUTION OF THE DIFFERENTIAL EQUATION 
the differential equation (Eq. 6a) is linear, solutions may be 
q by superposition of a particular integral and the solution of the corresponding — 
_-homogeneous equation. The form of the equation is such that ‘separable single _ 
series solutions of the Levy-type canbe obtained, 
- Solution m of the Homogeneous Equation.- —In view of the specified edge condi-— 


s (Eqs. 14a) a solution of Eq. 6a for a2= ve assumed in the form — 
= Only) cos =f Cc ne cos ) (15) 


= particular integral, the solution of Eq. 15 is genase ; general to — 
arbitrary boundary conditions along the twoedgesy=+b/2, . 
Introduction of Eq. ‘15 into > Eq. 16 for each oft thes series, 


We now on 


2 


; ; Substitution of Eqs. 17 and 18 into Eq. 16 gives the equation © — 


— 
| 
: 
a 
: 
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Lv 


which the quantities X1, My a and Up a are the functions of the 
shell parameters € andy: we 


| [v2 fan +4 
{1/2 


| [va «2. 


5 + Cg +Cge’ COS A: 9) 


sions, their form is such that a numerical tabulation may easily be made by 
7 means of an automatic digital ‘computer, ‘The use of numerical tables in d - 
4 To illustrate the use of Eqs. 21 for comparative shell sh shape studies, the 
i. 7 real parts 1 and X92 of the roots are plotted in Fig. 1 and Fig. 2 for various J 
DB ncqeer of the shell parameters € and y. Because of the exponential form of : 
2 the solution Eq. 15, 31 and X2 characterize the dampening of the edge disturb- 
a ances from the edges y = constant. Note the small values of Xo for positive | : 
values of y in Fig. 2, which indicates slow dampening of the edge mnnese 
for hyperbolic paraboloid shells, 
Particular Integrals of Eq. 6a: Surface Load. —Appropriate inte- 
grals can be cxpanting the load function in series. 


« 


— 7 _ Substitution of Eqs. 20 into the Eq. 17 gives the roots of the satin i 
— | 
— 
from whic obtained by substitution 
| 
| | 


ies 


7 


Tit 


Variable Surface Load Navier-Type Solution. — For various shapes o of the 


load function a particular tiss can be obtained in double Fourier series 
mn COs by cosvx. 


in which m and n are integers ai and and defined by 

The coefficients qmn are determined from the form of the load function, Sub-— 

stitution of the ‘Eqs. 24 into, 6a the coefficients 


\ 
“Ale.  & 
¢- 
j 
—_ 
a 1,—REAL ROOT x, OF THE CHARACTERISTIC EQUATION 
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FIG. —REAL ROOT OF THE EQUATION. 


f 
<4 
= 
4 a — 


we 
— 


_ supported edges along all four edges of the shell, It follows that Eqs. 24 are 
solution of the Ay problem given by the differential 


+ vot? By inspection it is seen from Eq. 26 that for Y > 0, that is is for hyperbolic 
‘paraboloidal shells, the factor (02, 62) may vanish for particular span- -cur- 


; direction, : a particular: integral can be obtained in single Fourier series form 


$n COS) x 


from which the stresses and ‘iain: are obtained by introduction into 
Eqs. 7. The complete solution of the differential Eq. 6a is then obtained by 7 
_ Superposition 1 of the solution of Eq. 23 and one of the particular integrals, Eqs. : 
24 or 


ary value problems, it is convenient 23 back into Eqs. 
_ express the statical quantities inthe shell in terms of the solution, In analogy 
_ With the methods used in the design of circular cylindrical shells, character- 7 


istic coefficients will be derived and a systematic solution procedure will be 


oan Characteristic Coefficients .—In view of the symmetry of the roots, Eqs. 21, 
Eq. 23 can be written in — ig form for each term of the series: 


— 


Ls ie... “Biegetheorie der Translationsflachen und ihre Anwendung im Hallenbau,” 
Hruban, Acta Technica 1953, Tomus VII, Fasciculi 3-4, 
> “On the Calculation of Shallow Shells” by S.A. Ambartsumyan, Prik, Mat. Mekh., 
XI, 1947, (NACA,TM 1425) 
“On the Calculation of Shallow Shells,” by WwW. ond D. A, Conrad, 
“Applied Mec Mechanics, ‘December, ‘1959, 


120 
from which the st ed by introduction into 
+ 
| 
ae _ the complete solution for simply supported shells reduces to the Navier plate — «ate 
a - solution, which shows that the hyperbolic paraboloid is an improper shell for a 
; 
: 
n 
— 
— 
29) 
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where the eight roots ‘roots from the the ‘Eqs. 21 are represented a by | two two princ 


5, 


(ating) 


- 4 and where C1- Cgare the eight complex constants of integration . By ‘definition 
the stress-displacement function is a real function, from which the following» 
"relations between the constants of integration results: 


a 


Introducing new constants, the so solution given in 1 Eq. 29 m: may be written in 
By inspection, it is seen from Eqs. 7 that most of the stresses and 
ments in the shell can be writteninthe form = 
= [H] Re| Cy fy « Hy +C3 Hy 1°" Hye (33) 
in whic] a [ H] i is a multiplier containing functions of the independent variable; 


x, and Ay, . and Ho are polynomials in the two principal ‘roots, Eqs. 30, which, 
by suitable choice of the multiplier [H] and change in the constants of integra-_ 
_ tion, can be made quantities of the order of magnitude 1. Using rite coeffi 
‘shell theory notation tl these naked will be denoted characteristic coeffi- 


= all the multipliers are reduced by a constant factor (1/2) Ehkge2, all of 
which represent changes in the arbitrary constants of integration ¢ only. Re- : 
_ duced values : are introduced for the roots of the characteristic. equation, that — 


nts 


(34) 


after numerical values have been introduced. 
Solution Procedure for Particular Boundary Value Problems. —The 


method for a boundary value problem given by the differential Eq. 6a, the bound- 


« 
o 
= Cy, Cg = Co, C7 = Cz, Cg = Cy (31) 
— 
| = 
a 
| 
| 
Characteristic coefficients Hj and Hg and the corresponding multipliers 
— : 
The characteristic coefficients for thedisplacement components uandvare 
= 


Fe 
"February 

ary conditions of Ea. and conditions the edges 


1. The actual edge conditions are stated along the edges y= = constant, 
_ The principal roots of the characteristic equation are computed from 


Eqs. 22 forthe actual parameter values and the characteristic coefficients a: are | 
the form formulas inTablel, 


"TABLE 1, —CHARACTERISTIC COEFFICIENTS 


2+ +18 4, 


4,7) 


3 


OR 


(eh, 


| 


integrals a are determined from Eqs. 24 or using the rel > rela~ 


4, The edge values of the homogeneous and the particular solutions are sub- a 
‘stituted into the edge conditions and the constants of integrationare determined, 
_ 5, Finally the stresses and the displacements inthe shell are found by su- - 
perposition of the homogeneous and the particular solutions evaluated for the | 


computed values of the constants of 


4 
ay 
- — 
— 


In dealing with edge | disturbances from two opposite edges, it often proves” 
‘convenient to consider each edge Time, | for the y= 
solution given in 32 32 reduces to 
ne for C2 | cos 


and the edge value of the wanes shell miles’ has the form _ 


typical sh 1 shell quantity is in the 


and the edge value at at ‘iy * 0 simplifies: to 


- [H] Re 


7 


the edge values are expressed in terms of the characteristic coefficients and 


the real constants Aj, By, Ag and Bg which in turn —— be determined from 


the edge conditions at y = - b/2, (y’ = =0). 

‘The value of the unknown statical « quantity at an arbitrary point in the shell 
is finally found by substituting the constants of integration back into Eq. 38. 

By expansion of the functions, Eq. 38 is transformed 

~ [a ag +ag fig +agf (41) 


cos (uy pry’) 


Hy{ - By 
| 


(36) 
&g 
i 
ual 
5 
(48c) 
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200037892 | 
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3 
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97964 
052750 
34795449 
34681079 
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Be 384376 


34228395 
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24830728 
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20385086 
24224882 
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(20099703 
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a 


244970060 
190962604 
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dl 
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Oe 
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— 

a _ TABLE 2,—ROOTS OF CHARACTERISTI 
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14880598 
le 847464 


16816884 
16788559 


14737708 
le 071478) 


1.693299 
16673124 
126654135 


14619304 

16603286 

16588100 
16573679 


1.538859 
16522574 
le 
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0.981725 
02941194 
00903834 
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00606810 
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00581767 


06535869. 
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"25 


quantities in the y- -direction, and he contributions must be subtracted from 
each other for the odd quantities, ee. 


‘DESIGN TABLES AND NUMERICAL EX EXAMPLES 
_— 2 foom of the solution “pene is such that it lends itself to tabulation as 
by Rudiger- listed influence coefficients circular 


— 
| 
0499 
— 
| 
1020 
— (21632 7 
1  Intreating the two opposite ed y, the effects from the two edge 
disturbances must be superimposed on each other. In superimposing such ef- 
4 
| 7 - load system on the edge is easily performed if the characteristic coefficients ‘ 
a are tabulated for various values of the shell parameters, It proves to be con- © : a: 
aa venient to choose the following unit line load system on the edge 7 “5 ‘a 

— My=0 My=0 j%My=0 jMy=1/[ 


3. —ROOTS OF CHARACTERISTIC EQUATION FOR 50 
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TABLE 3,- CON TINUED 


19756319, 
729254 
10704112, 
16680688 
10658805 
10619081 
16600991 
16583942 
1.567847 
19552624 
10524524 
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(00568496 
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06 494219 
(00481729. 
0.469881 | 
00458627 


00874081 
00849811 
00826419 
00521158 | Oo 803876 
00528456 | 00782154 
0542285 
0.546882 00721614 
00555302, 0. 702884 
00561564 | 00684833 
00567687 | 00667438 
06873684 06650672 
00634514 
06618940 


7 


00575511 


00513551 


04402097 


06562064 
Os 549101 


0.407806 

14455181 
(16445374 02428001 0.418078 
(14435982, 00418717 


where the symbol L. Pann that the stress resultant is reduced by its multi- 


hh an actual set of boundary conditions the internal aise a 


Wi 


that the shell and the edge member form a ‘compatible structure, 

—_ _ Tables for Parabolic Cylindrical Shells. —For y = Othe solution g given in — 
21 becomes identical tothe solution of the Donnell theory equationfor circular 
cylindrical shells as given by Ridiger-Urban, 5 provided ko is replaced by the 
inverse of the radius of t the circular shell. It follows t that the extensive tables | 
given by Rudiger- Urban can be directly adopted inthe design of shallow para-— 

bolic cylindrical shells. In applying these tables the tangential and transverse 
stress resultants and the flexural moment ts are directly replaced by the cor- 

: ‘responding quantities of the circular shell. The displacement component w w in 
the z-direction is replaced by the radial displacement vz,z, and the first order 
_ derivatives of w with respect | to x and y are replaced by the | corresponding — 
slopes for the circular shell, Also note the quantities 
are all reduced by 1/Ve by Riidiger-Urban.® 

“OM Tables for Hyperbolic and Elliptic Paraboloidal Shells .— 


00524549 


“a and elliptic paraboloidal shells by assigning various values to the shape a 
rameter y between + 1 and - 1, A general program for such table computation | 
by means of an automatic ‘digitalcomputer isin progress. An assembly of such © 
tables for a suitable choice of parameter values will be the subject of a sep- : 
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— 0087 
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(1614 
1e17 
1629 
(45) 
_ 
— a 
&g 


oF 


0002 | 76823723. 
~ 0903 
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49254869 
34993643, 
- 
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39440416 
31304967 
| 36185870 
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(24822766 
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(16704601 
10676708 
1665097? 
10627146 
10605021 
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EQUATION FOR 
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Be 
34182833 
30035942 
22905883 
20789610 | 
20589635 
2e 502691 
20422816 
20349073 
20280693 
0157559 
2101812 
20049403 
22000000 
10953314 
15909093 
16867119 
827198 
10789160 
10718146 
10653057 
10622470 
10593069 
10564775 
10537516 
“test 
10485846. 
te 461321 
10437602. 
10414642 
__10392400 


6b 


19253856 

(10109442 
10067654 

06 991613 
“00956869 


00771469 | 00924041, 


— 
| 
— 
— @ 
q 
im 
2 
avd 0057 | 
— 


TRANSLATIONAL SHELLS 


210530286 | 
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TABLE 4,-CONTINUED 
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‘Tables 2, 3, 4, and 6 give the ‘roots of the equation for two 


“hyperbolic paraboloids, cylindrical shells, and two elliptic paraboloids, that 


‘For t the case of a rotational paraboloid, that is for y = - 1, the differential 


Eq. 6a reduces to | to the weel-| known form ang 


the solution degenerates to a sum of the solutions for the free 
for the plate on elastic foundation, The procedure previously outlined d does not 
apply directly to this case, in which the biharmonic operator vanishes. The 7 
; multiple roots of the biharmonic equation are represented by the first princi- 
- pal root, (Eqs. 24) where by b becomes identically equal to zero. Table 6 is in- _ 
_ cluded for illustrative reasons, The roots dt and Lo represent solutions for 
“the rectangular plate onelastic foundation, 
Numerical Examples. .—The ‘solutions obtained | are particularly useful for 
‘comparative studies of shell shapes. - To investigate the effects of change in 
shape, four different shells have been analyzed, two hyperbolic paraboloids, — 
one cylindrical shell, and one elliptic paraboloid, which had the following val- a 


ues of the shape parameter y vy + 1, + 1/2, 0 and - 1/2, respectively. The pa- 


Jf 
0684 
ii 
00S 
1023 
163 
| 
§ 
| 
2 


0201 
0.02 


16401803 
10412357 _ 


16414035 
16414193 
10414183 

106414181 

10414145 
16414116 

14614076 | 

10414024 

16413959. 


16413780 | 


16413663 


19413522, 
16413164 


16412392 
164122061 

411689 


1410296 
6409731 
0409112 

408837, 
0407704 


16405137 
10406156 
16403110 


1» 398285 
15396920 
16394007 


16382080 
ts3T6286 
16370186 
(19363854 
10357357, 


| 


0.000556. 
06001270 
0.001729 
00002260 
0.003533 
0.004273 
0.005085 
04005966. 


0.009023 


Os 5016933 
02018471 
00020066 


00026983 
00028832 
02030723 
00032652 
0.034616 
00038633 
0040678 
00042741 
02044819 
O06 0046907 


00053190 
00055277 
0.057353 
059416 


0067456 


0 073182 


00083602 

00092451 

04096269 


0.099693 


~ 2666s 


700710802 | 


356355379 


U 
174677699 | 
146142184 
112785191 
104101649 
84839018 | 

| 
Sa 071429 

893184 

54440085 
_5e 051769 


420936 
= 161281 


34930554 


724197 


3.538563 


(UTUT— 
39218199 
30079063 
20951633 
20834513 
20726525 


24534069 
20447999 | 
20367810 

20292945 
24722912 

28t 
24095676 
29037753 

16983219 | 
10931805 
de 

10794010 
10752912 
10676987 
16441883 
US™ 
10439321 

es 
16311314 
19258917 
10712768 

16171995. 

16135872 


700710555 

170677688 
140142176 

100101645 
Be _ 80899013 


428702 
50893149 
50440032 
56 051693 
46420787 
40161081 
30930289 
39723850 
30538115 


e U 
30217483 
30078171 
950534 
20833171 
20724900 


20531741 
30445262, 
20364568 
20289154 
| 20218506 
20089820 
20031057 
10975597 
| 10923169 
10873531 
de 781771 
10739274 
10623424 


10491435 


0260938 
10198826 
4 0142224 
10090344 
10042545 


— _ TABLE 5,—ROOTS OF CHARACTERISTIC EQUATION FO a 
| 
| 06 10 | 
0018 | 06011400 — 
| 
— | | 
| 
8830 
0031 
— 

| 0969 

| — 


lel? 


1623, 


TABLE 


19324305 
16317853 
19321515 
16305303 

U 


14293303 |, 
16281904 
(19276436 


le 260953 
de 256092 


14251375 


19246799 © 
16242361, 


14233879 | 
16229827 


le 225896 


“TRANSLATIONAL SHELLS 


00105435 
00107798 


06115499 
117078 
‘Be 117623 
00118027 


T8308 


00118466 
00118526 

00118494 
06128380 
06118193 


00117628 


(06116853 
gia 


> 


NTINUED 


14075247 


10049800 

10006789 
| 


06957867 
(00944866 
00922850 
0. 913575 
04897989 | 
06 891415 
00885628 
06 880517 
00876019 
TT 
(02865660 


06863095 


00957182. 
00918827, 
00849255 
817569 
00759475, 
00707465 
06 
00660638 


00618289 
00598614 
00579858 
00561965 
00544885 
04512988 
00483839 


rameter € was chosen € = 0,30 for all. shells. All shells were square in plan | a 
i =b), and all shells had free edges without edge beams along the two —e 


All shells had load lof, = 
ib. For this load distribution a particular integral is found from Eqs. 27 and 28, 


which by substitution into Eqs. 7 gives 
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Nyy = = Nyx 
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€ cos Ax ) 


were determined | 


from | Eq. 39 and ¥ were e corrected for effects ene the opposite edge by means — 


Eq. 41, characteristic coefficients were from 
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ratio was chosen as v= 0.25, 
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to, and ") from Tables : 2, 3, 4, and 5. The con- 


- stants of integration were obtained from the boundary conditions of Eq. 47 by 
_ introduction of the particular integral of Eqs. 48 and the edge values from the 

For illustration the transverse distribution of the flexural moment Myy and 
g the longitudinal normal stress Nxx resulting from the first term of the series _ 
_ expansion of the load function are shown in Figs. 3 and 4, respectively, = = 
The numerical computations involved will be shown for the hyperbolic para- 
- ‘poloidal shell with y = 0.50. Only the first term of the series expansion will i 

‘The roots of the characteristic ‘equation are from Table 3 for the 


= = 74 
The characteristic coefficients are then computed by means of the formulas b * 


"listed in Table 1 and are given for a number of stress resultants in Table 7, 
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‘The distribution functions f1--f4 are determined from Eqs. 42 anda are given» 


that effects from the edge disturbance at the ineiie edge must be considered. 
The edge values at y = = Oresulting from from the | edge disturbance at -_ b are found 
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-1,13595 +i1,54917 .725278- 128372 


ll .034552- -i3, 519544 0. .584552- io, 186210 


13595 A, - 1.54917B, + 854578 


— 
&g 
9 = -0,3636 A> + 0.07266 Bo | mas 
= -0,07044 Ap- 0.2852 Bo} 
_ Since the boundary conditions of Eqs. 47 are symmetric, the total edge val- 2 
oe. : ™ at vy = 0 are obtained by adding Eqs. 49 to the values obtained from Table — a 
COEFFICIENTS FOR = 0,50, €= 0.30 

Nyy 
4 
— 
‘2, Introduction ¢ 
boundary conditi — 
a ... 60) 


— 
solution of which 


= -1,02956 q’ 


The s stresses and are then obtained by of the con- 
stants of integration into Eq. 41. The stress resultant _Nxx and 1 the flexural 
NXX 


Nyx = [ex] 7.630 fy + 4,996 1.755 fg + + .628 8 


[Myy (1.56 563 f, + 2,311 fp - 0.188 fg + 0.604 f4)q’ 


distribution in the y- -direction is most in tabular 


as shown in Tables 9 and 
TABLE 8. _—DISTRIBUTION FUNCTIONS Ys = 0. €= 0. 30; a=b 


f 


3 
4 
6 
8 
9 


af 
Analogous computations must be made for n = 3, 5--. . Note that the ‘shell 
h?3—ee € is a function of n. For uniform load the series converges rapidly 
and a good approximation is obtained by the first term only. ee 


CONCLUSIONS 
aa Levy- -type solutions of the aio of bending of thin shallowtranslational _ 


‘shells are obtained. The solutions are expressed in termsof two parameters, “ 
one shape parameter y = or y= in which cj and are 
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application in actual boundary value problems is 
of tables of the roots of the characteristic equation are included, and the use 
of tables in application is described. Graphs showing stress distributions in a 
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shells that have method ome is 
- shown that the dampening of edge disturbances is slower in hyperbolic para-_ 
_ boloidal shells than in cylindrical shells and elliptic paraboloidal shells, which 
general indicates that edge disturbances produce relatively more significant 
bending ¢ in n the hyperbolic paraboloidal shells. 
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alate. —After the final copy of the present paper was submitted t to ASCE for 


_ publication, it came to the writer’s attention that the work presented is paral- — 
in part, to developments made by A, L, Bouma. 
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WIND- )-INDUCED VII VIBRATIONS IN IN ANTENNA MEMBERS 


7 By William Weaver, co M. ASCE © 


This paper the of wind- tunnel tests of wind- 
tions in the circular aluminum tubular members of large radar antenna space 
‘frames. A curve for the RMS von Karman lift coefficient is established, and a 
pp gpeetectorede means for evaluating the effect of self-amplification of vibra- 
tions is presented. A system of ‘spoilers is also investigated and 


recommended for applicable cases, 
— 


and Background of the Pr oblem. —A steady flowa ona 
-lindrical structural member induces lateral vibrations in the member due to 
the formation of vortices on alternating sides of the member, — These vortices 

are’ formedand shed from opposite sides of the member ata a regular frequency, 

; depending on the velocity of the wind andas a result alternating lateral forces” 
are exerted on the member by the fluid motion, When the frequency of vortex — : 
_ shedding (and hence, the frequency of the alternating s forces) is approximately 7 

a so to the frequency of one of the normal modes of vibration of the member, _ 
a condition of resonance results, If the structural damping and the stiffness of 

‘member are small and the wind large amplitudes of 


"may be. considered to be self-exciting for alimited range of wind velocities sat 


_ Note,—Discussion open until July 1, 1961, To extend the closing date one month, a 
_ written request must be filed with the Executive Secretary, ASCE, This paper is part 
_ of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 87, No. EM 1, February, 1961, 
1 Asst. ‘Prof, of Civ, Engr., Stanford Univ., ‘Stanford, vou 
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and around the resonant (or “critical” ) wind velocity. This 


duced by the wind are not catastrophic in nature, Oscillation persists for a. 
‘certain of wind ‘Speeda above the critical and a ‘maximum 


the stiffness of | the so that its 


vette? is above the range of moderate winds; 


_ 2, By reducing the effective length of the member through the introduction | 


_ 3. By use of damping devices to restrict the amplitude of vibration; or 
4, By cores “spoilers ” to the member that serve to i ly the flow 


destroying the cause of the vibrations, 

In the past, the first three methods of suppressing wind-induced vibrations have 
been used repeatedly. — The solution to this problem by the method of applying» 
spoilers to the surface of the member has not, in general, been attempted in 
civil engineering structures because of a lack iene concerning wal 
‘design and their effectiveness. 


This _ investigation originated v with the problem em of wind- -induced vibrations 


have been attributed to the action of von Karman vortex forces. The failures” 
occurred in slender tubular members subjected to moderate winds, and they — 
were observed to be fatigue- '-type in weld- heated areas at points of 
Notation.— The letter symbols a adopted for use in this pa paper are defined where 
they first appear, in the illustrations or in the text, and are arranged alpha- 
for convenience of reference, in Appendix I, 
Object.—One of the objectives in the current investigation is to obtain an 
: experimental curve ar! the variation of the von Karman coefficient, Cy, with» 


imately from 6. continuous curve for CK in the 1 region of interest 

is a basic necessity for the establishment of design criteria for wind- induced | 

aia A second objective is to study the self- amplification aspect that concer 

_ izes wind-induced vibrations of this type. With sufficient experimental data, 

involving variations in all of the parameters affecting the response, a —_— 


empirical solution to the problem can be applied, 

_ Strouhal numbers, §, are evaluated in this investigation for the purpose of 

comparing with well-established values previously published in in engineering 

- Another purpose for these experiments is to test and optimize a system of 

_ spoilers attached to the surfaces of a variety of member sizes. The function 
4 of these spoilers is to interfere with the periodic circulation and formation of 
vortices around the periphery of a circular cylinder. The measure of their 
- effectiveness, Se, is the ratio of the von Karman coefficient | fort the cylinder 

_ with spoilers, Con to the von Karman coefficients for the bare member, CB: - 
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_ — This investigation is limited to a study of tubular aluminum mem- 
“eee in the form of nenaaabtaneennintal The types of test members considered © 


are as” follows; 


‘Stationary rigid- 
Spring-supported rigid- body pn 
End- ‘supported flexural members - fixed; 
4 End-supported flexural members - partially fixed; 
5. End-supported flexural members - pinned; and t 


The spoiler configuration selected for these consists of 
diameter tubing attached to the surface of the member and wound in a helical ; 
_ pattern over all or part of the length of the member. This particular config- 
uration was chosen because of the ease with which windings can be installed 
_ on a cylindrical member and because of the fact that a helical patternof spoil- ‘ 


ers is not sensitive to the direction of the wind, The variable uneamaenes for 


_ optimization of spoilers of this kind were the followings a ote 


a. Number of windings; 
Pitch of windings; and 


an: - Portion of length over which windings are 
_ The appropriate range of R was obtained by means s of ntieeine standard- 


diameter test members and cylinders | to the full range of wind speeds in a wind 
tunnel having a velocity range of 10 mph to 100 mph, 
a In tests on vibrating members, stiffnesses and structural damping were var- 
ied sufficiently to allow the effects of these parameters on the the response t to be 
_ The test program was divided into three main phases: Phase A ‘consisted of 
tests on stationary cylinders mounted in inthe wind tunnel with a sensitive trans- 
ducer built into the supports | insuch a manner that direct measurements of the 
reactions due to the lift forces acting on the cylinders were obtained. Spoiler 
: optimization was carried out on the stationary cy linders and confirmed in the 
7 later phases. Phase B dealt with the response of spring- -supported cylinders 
oscillating as rigid bodies in the wind stream, Phase C was a study of the re- 
‘sponse of flexural members vibrating in bending modes under the influence of 
a steady wind, Each of these phases contributed information on the lift poles 
r 


3 ficient, the self- -amplification phenomenon, the Strouhal number, and the he spoile F 


RESPONSE CONSIDERATIONS 

Response of an Elastic System to Vortex Forces (Ww “ithout Self- -Amplifica- 
Sinusoidal Variable Forcing Function,—The commonly expression fc for 
the alternating lateral force due to vortex is (1):2 
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which’ Fx is 3 the von Karman denotes the stagnation is 
the projected area, f y represents the frequency of shedding vortex pairs, and 7 
tis time. For a given velocity this idealized forcing function varies a 
ally with time and has an unvarying and a discrete 


in which y is the | amplitude tT { vibration ft an elastic system, d represents the 

diameter of the cylinder or member, C; isa constant for a given system, 2 
Gar the ratio of frequency of forcing function to natural frequency, k is the 
e spring constant for an elastic system, and 6 represents the the logarithmic damp- 


ing decrement, or rate of decay. 
_ Eq. 3 may also be applied to the case of a flexural member v vibrating in its _ 


‘fundamental mode, in which case the spring constant, k,,, may be approximated 
as the total load acting on the member divided by the maximum deflection of 
the member when that load is } applied re 


, Random Forcing Function.— For R above approximately 3 x 102 the ampli- 

' a of the lift force is not constant for a given stream velocity but woregneadl 
108 increasingly random as R increases, At an R-value above approximately 3 x 
105 eng vortex- shedding frequency also becomes random as the Separation point — 
‘moves rearward on the cylinder. Therefore, for the purposes | of ‘computing 
yy the response of an elastic system, there are two important types of forcing — 
functions to be considered: (1) 3x 102 <R< 3x 105, sinusoidal forcing func- 
tion with discrete frequency but random amplitude; and (2) 3 x 105 <R, har- 
- monic forcing function with random frequency and random amplitude, 7 
Bey _ Response at Subcritical Reynolds Numbers.—The range of R below approxi- - 
‘mately, 3x 105 is of primary interest in this investigation because most prac- 
_tical- sized antenna members subjected to moderate wind velocities fall into 
a ph In this case the force e amplitude i is random with a ‘a steady mean v: value of zero. 
- The response of the system depends on the energy input from the vortex forces 
averaged over a periodof time. Now, the energy input per cycle is proportional 


7 to the product of the force and this ‘amplitude of vibration. But the amplitude | 
_ itself varies directly with the applied force, so the energy input per cycle var- 


ies as the square of the force: 


Energy input per cycle «Frye FR atch 


/§ ‘Therefore, the mean square of the force determines the energy input, and 
the response must be computed on the basis of the standard deviation from the 
mean, or the root-mean-square (RMS) of all of the random deviations, We- 
may define the non-dimensional lift coefficient, Ck, as as the RMS- vee of CK 


p CK sin 2 
a Ap esin27f,t 


and compute the ‘response from | 
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WIND INDUCED VIBRATIONS 
Response at Supercritical Reynolds: Numbers.—Supercritical R- -values are 
_ of secondary interest in this investigation but are of primary interest for struc- 


In this case the force amplitude has a mean value of zero, and both the am- 
_ plitude and frequency are random, We must, therefore, deal with power input 


rather | than energy input per cycle, because the poe depends on the frequency 
of the excitation as well as the force amplitude: 
Power o y ty « f 


> _ resulting plot of power versus frequency is termed the power spectrum, and : 
the ordinate associated with a particular referred to as the “power 
= ‘spectral density” (PSD) for that frequency (2)(3). A frequency of high spec- 


= tral density will produce a high response at resonance, a 768 


For the stochastic process under consideration, if. the power spectrum “te 
expressed as a function of 


Om 


then the mean square amplitude of the forcing function may be expressed in 


terms of the mean square lift coefficient and the integral of the power spectral 


‘The mean square value o of the response may then be computed from 


. Under the as assumption that there are no : no phase « differences along the length 


of the cylinder, and not considering the effect of self-amplification, Eq. 11 may _ 
_be evaluated | with respect to the power spectrum in a narrow band width at a 
frequency corresponding to the natural frequency of the elastic system (4), . 

_ The mean in square response becomes (approximately) 


o(S). wee 
“High spectral densities at low $-values would indicate the probability of 
dangerously large oscillations at fairly high wind velocities, ve Trew 


Self-Amplification of Vibrations. = 
Incremental Load Factors.—When an elastically supported cylinder or a 

_ flexural member is vibrating in resonance with periodic vortex forces, the vi- 
bration itself causes an increase in the magnitude of the lift. The additional — 
aerodynamic fo force is probably associated with a periodic ae of the > points 


— 

— 
— 
a ; The power input must be evaluated interms of the whole spectrum of frequencies © : | = 
____ involved, In order to do this, some characteristic function for the power with a 
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time, the vo vortices created in the wake are age and stronger than inthe « case 


ae B. Steinman (5), citing the results of experiments at the David Taylor 


4 - Model Basin, proposed the fe etn that if eddies are shed at or near each 


(13) 


_ in which h isthe width stieiniians the centers of two rows of eddies behind a non- 


oscillating circular cylinder 
sii For ‘stability, th the distances between vortices in each row must also increase 


in the same ratio. It isassumed that the drag, the: circulation, and the lift must 
all increase proportionally, At the same time, the eddy frequency and $ are 


‘There is an upper limit of stream velocity, vp, beyond which this self-con- 


trolling action cannot occur, At this point the amplitude drops abruptly to a 
— lower value determined by the en factor or corresponding to to the fre- 


s 


The ratio of vp to Vor (the critical wind velocity, in which fy =f, and fp, i 
the natural undamped frequency of an elastic system) m must be the same as the 
ratio of dimensional increase in the vortex street: ‘a: 


For the case of zero o damping, Steinman | ane an upper limit of velocity 


Vp © 1.45 


7 te ratio of the self- amplified lift force per unit length of an po: A 
7 cylinder or flexural member to the lift force per unit length, | if there were no 


amplification, ma may be e defined as the | “incremental load factor” (ILF): 


ancecillating cylinder th the ILF- value i is equal to to For a a 2 vibrating 


vibrating under the of a wind id velocity somewhat than the crit- 
: ical velocity, the aerodynamic forces exceed those acting at the critical veloc- 

‘ity for two reasons, First, the forces increase in proportion to the velocity 
_ pressure h head increase (v/ver)?. - Second, the forces increase a to the self-— 
pence for both of these factors we may define the “ incremental deflection 
factor” (IDF). The IDF-value is the ratio of the peak amplitude to the hypo-- 

_ thetical amplitude at the critical velocity. The latter amplitude must be rec- 

; ~ ognized as hypothetical because resonant vibrations of large amplitude cannot 
be attained the for a cylinder equals: the nat- 
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- frequency | of the body. In this instance, as soon as oscillation begins, the 
ects frequency decreases and goes out of phase from the natural frequency, 
Fora cylinder oscillating as a rigid body the ILF- value and the IDF- value 
are determined by the damping, the spring stiffness, and the lift coefficient. _ 
_ a vibrating flexural member, however, the ILF-value varies along the | 
= of the member; therefore, the IDF- value is a function of the shape of — | 
‘the dynamic deflection curve and the phase relationships of the forces acting 
soi phase differences along the lengths of flexural members are ignored, the 7 
deflections of all vibrating members may be approximated by fourth-degree _ 
_ parabolas. +The incremental loading on a member may then be assumed to vary 
as a fourth- degree parabola, and the incremental deflection factors for vari- 
ous types of members are then found in the following manner: 
For a simply supported d member vibrating un under the ‘influence of the wind 
elocity that causes peak ge 


384 Ta T 


in which Ymp is the amplitude of vibration at point of maximum deflection on 


_a flexural members when v = = Vp, including the effect of self-amplification; we 
denotes the effective lateral hows per unit length due to vortex shedding on = 
a cylinder with v= Vp» ] is E the modulus of elas- 


ve 


5 wt 


3d 


ve) =n 1.779 x 10- ° Cx Vv 


22 into Eq. 


lowing which A is the cross-sectional area): 


Pinned Members: 


— 
a 
4 
| 
4 
hn 
= ee eeeee : ( 1 8) 
definition 
_in which w is ty of air. Also, 
— = 96,400-5 + (22) 
in which Is the radius gyration, Substituting Fqs. 
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_ Aeroelastic Suppression Factors.— — The ID IDF-values in Eqs, 23 | through 2 26 }de- 
pend, among other things, on the ILF-values. The ILF-value depends, in turn, 
on the following three elements: (1) the damping, (2) the stiffness, and (3)the _ 
relative magnitude of the wind forces, a These parameters may be be “lumped” in in- 
toa dimensionless number that will be called the “aeroelastic suppression fac- : 
7 tor” (ASF). The ASF-value is a reciprocal measure of the relative ability pot 


vibrating systems to amplify their own vibrations, > 


Expressing ig Ke al and Fx in terms of m, $, and —— 

ve 

in which m is the mass per unit length a spring- 
flexural member, The ASF-value for a flexural member vibrating in its fun- # 
damental mode may be expressed in a manner ‘similar to that for a spring- 
supported cylinder. The stiffness, k,,, of the flexural member, however, will 
be defined as the total uniformly distributed static load that produces a maxi- 

_ mum deflection of unity. Accounting for the fundamental mode shapes and oil 


a _-Tesponding member stiffnesses in this manner, the following ASF- values: are ‘a 
m8 28 
ASF = 5,28 x 10° - 
5.44 x eee ee eee eee, 


WwW 


— 
| 
q 
> 
f 
f 
— 
— 


“WIND INDUCED VIBRATION: 


i Values for ILF and ASF may be obtained experimentally, and a plot of ILF 


_ versus” ASF can be constructed, The resulting curve serves as a useful tool | 
for the evaluation of Eqs. 23 through 7 and hence, the degree of self- ampli- 


_ Experimental The experiments inthis investigation were con- 
; ducted in the low- speed wind tunnel in the Aeroelastic Laboratory of the Mas- 


sachusetts- Institute Of Technology, Cambridge, ‘Mass, This tunnel is 3 of the 


clo osed-circuit _ having abd ft by 7 } ft rectangular test section and a ve- 
lo ocity | range of 10 mph to 100 mph, . Screens in front of the transition section 
served to reduce the scale of turbulence at the test section to approximately © 
5% Tunnel controls» were and blockage effects (of test members 


determined using a 


an alcohol manometer, = 


A Sanborn oscillograph was used in the experiments whenever feasible, In 
“cases ir in which the frequencies of the electrical signals exceeded the capacity | 


_ of the Sanborn, a Heiland oscillograph was used, 


‘Tests onStationary Cylinders (Phase A).—Fig. 1 shows the base layout and 
- the ania transducer connection for the Phase A tests (with one of the _—— 


4 


4 
1.—PHASE A, BASE LAYOUT (3-IN, DIAMETER CYLINDER) 
4 
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Calibration---Sanborn Oscillo raph 


(10 lb ea. 


way at x20 attenuation = 20 mm amplitude. 


A 


= 20. 0 


A 


Example Sanborn Oscillograph Record for Stationar 
8 Windings Paper speeds 10 mm per sec.; v = 


Atten. x4; Fx = 2, 4 1b; Cx = 0.519; 18.0; 


= 30.0 mph; 


0. 184; 


; 


Example Sanborn Oscillograph Record for Stationary with 
$poilers Paper ‘Speed 100 mm per sec.; 30. O mph; — x45 
0.15 1b; Cy = 0.035; R= 
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INDUCED VIBRATIONS 


guide plates removed). A “dummy” connection of similar construction was 

_ provided at the upper end of the cylinder so that the force reactions at the top “ 

and bottom connections would be approximately equal, 
The force transducer used is basically a compression force load cell with a 

25- -lb working range. . The load cell was modified for measuring alternating 

_ tension and compression forces, however, by precompressing ang force-re- 


range of the instrument was + 12.5 Ib, instrument was cal- 
ibrated with static loads previous to mounting it in the tunnel, and it was again 
checked after mounting by means of applying static lateral force to the cyl- 
_ The bare cylinders were mounted in the tunnel and subjected to steady wind 
speeds in 5 ‘At each velocity an oscillograph record was taken, 


| 


‘ 
g 
— 
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an example | of which is shown in Fig, 2(b). For | certain of the bare cylinders ~ 
RMS galvenometer voltages were recordedusing a voltmeter-analyzer arrange- 
ment, | This set-up consistedof an AC electronic voltmeter in combination with 
a DC voltmeter shunted across the recording galvenometer. The former meter : 
_ was used to determine values of Cx, and the latter meter provided a means of 7 7 - 
load calibration, The voltmeter readings \ were taken with 


ay. potions the course of spoiler optimization the fact was discovered ‘that -con- 
figurations of many small-diameter windings had the effect of increasing the 

Es ‘lift forces ‘on stationary cylinders instead of suppressing them. Moreover, the 
«lift forces” not only larger in ‘amplitude, but they were more regular in 
frequency and also more . uniform in amplitude than for the case of bare cylin- 7 
ders, Furthermore, this force regularity transcended the critical Reynolds 

_ numbers for bare cylinders Fig. 2(c) demonstrates the e effect of many small 

windings applied to a stationary cylinder, 

Peete on Spring-Supported Cylinders (Phase B).— The same cylinders were | 

_ used in Phase B that were used in Phase A. The Support system was trans- 
formed by removing the lateral restrainst top and bottom so that the cylinders 
could oscillate freely on their spring supports, Twinaccelerometers were at-_ 
tached at both ends of the test cylinder to complete the Phase B set-up, Vari- 
able spring constants for the supports were obtained by using different lengths 
ss thicknesses of the cantilever support strips. Bert strips were clamped 


tensioned with a torque wrench, Fig. 3 shows the 4- in. ‘diameter spring-: 
cylinder of Phase B with optimum spoilers attached, = 
‘For each test on each bare cylinder accelerometer readings were recorded 
- for a series of wind speeds, ‘and the resonant wind velocity was carefully a ad- 
_ justedin the critical range sothat the peak amplitude of vibration was observed 
: in each case. Fig. 4(a) represents an example oscillograph trace for a bare 
cylinder oscillating at maximum amplitude (note the steadiness of the vibra~ — 


Previously Spoiler configurations ‘dramatically re reduced the re- 


‘cillogragh record of Fig. 4(b) with that of Fig. 4(a)(the effect is even — 
- than a first glance inc indicates | because different gains were used for the the two rec- — 


a ‘Tests on Flexural Members (Phase C).—The types of members tested v were 


Fixed in. to 3 in. 
‘Partly Fixed | in, to 3 in, 
a Pinned 1} in, to 4 in. 
Canti Cantilever in, to 3 in 


Cantilever members with 1/d = 30 and 1/a = 20 were wneeia by. cutting off” 
one end of the fixed members after the M1 tests had been campleted. In this. 
manner a i total of s sixteen different mem mbers were tested using only ten shop _ 
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ee Example Sanborn Oscillograph Record for 6"$ Spring-Su 
Oscillating at Maximum Amplitude rf 
Paper speed 10 mm per (Sec.; Atten. = 0. 022; 


29.8 mph; accel, = 00g; y/a 0. 036. 


(6) Example Sanborn Oscillograph Record for 6"¢ Spring-Su 


with Spoilers 4-3/8" at 12d, Oscillating at Maximum Amplitude 
Paper speed 10 mm per sec.; Atten. zi; 0. = 533 


= 15.05 cps; v = 31.3 mph; accel. 07g; = 0.0005. 


Example Sanborn Oscillograph Record for bl | Vibrating at Max. 
Paper speeds 10 and 100 mm per sec.; Atten. x 20; 


5 = 0.013; fa = 20.0 -_ vp = 13.21 mph; strain = 300 p in/in; 
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REPORTED EXPERIMENTAL VALUES OF Ck 
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FIG, 5,—PHASE C, END CONNECTION FOR 2-IN, DIAMETER MEMBER 
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‘Fig. 5 shows: an example end ‘connection : for Mlb as it was mounted hori- 
_ zontally in the wind tunnel, The bolted connections were tightened uniformly 
_ with a torque wrench, and the guide plates (one of which is removed) were cut 
to clear the members by ¢ in. Bonded strain gages (SR- -4) were used in the 
tests of Phase C to measure vibrational fiber strains, 2 — 
Bare flexural members were subjected to the range of wind velocities as in 
the previous tests, and fundamental modes of vibration were readily ‘developed a 
(Fig. 4(c)) at or near the predicted critical wind speeds. In the more flexible _ 
- members the secondand even the third modes of vibration were developed with © 
similar ease, The layouts of strain 1gages were oriented to the first mode, how- 
ever, so no quantitative results are reported for these higher modes. 
_ Optimum spoilers were applied to the members over various portions of — 
their lengths, The length increments tested 100%, 


20% of the lengths of the 
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Results. 
Coefficient.— 
: tion are presented in 1 Fig. 6, together with values of Cx reported by previous _ 
_ investigators. Table 1 summarizes published values for Cx whether 
“they are experimental or analytical, 
_ The three curves for standard CK from each of the three phases merge into” 
each other fairly well except where the curve for vibrating cylinders descends 
x 105, Moreover, the Cx-curves show 
“fair correlation with curves by Fung and Fujino for RMS Cx shown in Fig, 6. | 
“4 The upper Cx-curves for cylinders fitted with numerous small-diameter 
_ windings delineate the high lift forces observed in these cases, The peak CCK 
; curves for eight and sixteen windings correspond closely to the values given 7 
so Bd McGregor and Goldman and the curves given by Fung and Macovsky for ex- 


treme values of CK It appears, then, that applying a pattern of ‘small windings 
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WIND INDUCED VIBRATIONS» 
‘tot the surfaces of cylinders inthe present investigation has the effect of raising 
the peak values of Cx to approximately the same values reported by several | 
— A possible explanation for this phenomenon is that for some > 
reason, spanwise correlation of vortex shedding is improved by this type of 

surface roughening. | ‘The peak ‘Cx-curve for bare cylinders probably depends 
_ to some extent on the turbulence characteristics of the wind tunnel for its rel- 
Ss location, Neither this curve nor the other peak eak Cx-curves h have any prac- a 

tical value for the purpo purpose of computing the response of a vibrating member, 
It is the that assumes great importance inthis respect. = 

_ Self-Amplification.—The ILF-values computed for the cylinders and the 

_ members at resonance are plottedagainst ASF in Figs. ‘T and 8. Similar curves 
Ss Vp/Ver are practically identical with the ILF-curves, Thus, the « supposition — 
: that ILF = Vp/Ver is verified, and the ASF-value proves itself to be a valid 
_ basis for determining the degree of self- -amplification. The ILF-data for the © 
7 _ vibrating members is more scattered, and the curve is definitely lower than 

: ‘for the cylinder tests. The latter discrepancy is probably dueto phase areca 


~ ences along the lengths of the members that aa not arise along the lengths o 


cylinders oscillating as rigid bodies. 
_ Strouhal Numbers,—Curves for $ obtained in this investigation are come 


: pared in Fig. 9 with a curve obtained (15) by B. Etkin, et al (shown dashed) 

that represents a summary of the works of several previous investigators. _ 
general, the curves check farily well, and the maximum difference in ordinates 
_ is approximately 5% 
Spoiler Optimization. — Fig. 10 demonstrates, in a general manner, the in- 

3 ‘fluence of the number of windings on spoiler effectiveness. Configurations with | 
four windings repeatedly gave maximum suppression of lift forces. 2 ae 
«Fig. 11 illustrates the influence of winding diameter on spoiler effectiveness” 
for the 3-in. cylinder. The knee of this _ of curve occurs consistently at a 


= 1 for small-diameter windings and appear to rise only slowly for large —_ 
The relative insensitivity of the phenomenon to spoiler pitchis sue in Fig. 
12 forthe 6-in, cylinder. A wide effective range exists froma pitchof approxi- 
_mately 8 dto about 16 d. The most effective pitch appears to be at about 12 d. . 
Thus, optimum spoilers to be 4 windings of size or 
_ Sample response curves for | a bare test member and for the member with 
spoilers over varying portions of its length are presented in Fig. 13. For most 
of the members tested, vibration amplitudes were suppressed to approximately 
10% by optimum spoilers wound over portions as short as 20% of the member | 
lengths. Aspoiler length of 0.40 gave the best results in most cases. Spoil-— 


ers were most effective | for relatively stiff (fixed end) members at R-values — 


in the vicinity of 105, 


Correlation. —The question spanwise correlation is not consid- 
“ered in this paper. D. M. Macovsky (3) produced visual evidence that the flow 
picture at subcritical Reynolds numbers is not two-dimensional in its nature, 

_ The extent to which oscillation of a cylinder aids spanwise correlation is an — 


question that warrents further investigation, 
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12, INFLUENCE OF PITCH ON SPOILER EFFECT TIVENESS 
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‘The ability of vortices to conform to the m 
_ demonstrated by the ease with which higher modes of vibration are iene, 
_ Inthese higher modes, vortices must be shed simultaneously from opposite 
a sides of the member along different portions of the span, and yet the flow pic- ; 
ture readily adjusts itself to accomplish this complex feat. 
_ Cautions on the Use of Spoilers.—Spoilers are not recommended for Reyn- a 
olds numbers below 1 x 104or above 4x 105, because it was in this range that 
7 spoiler optimization was carried out. - Fluid flow has greater ability to over-_ 
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“ride protuberances at “the e lower R- At high the lift 
‘forces on bare (cylinders: are quite erratic without the > addition of spoilers, and 
Pha: 
critical R-values with and without spoilers, 
pay Spoilers are not recommended for members with diameters s smaller t than 
2 in, or for cantilevers with Va greater than 20, This statement is based on 
the fact that test members having small diameters and large nondimensional - 


= amplitudes o of v vibration yn during bare runs demonstrated an ability to persist in a 
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CONC LUSIONS 


from Tests.—The experiments conducted in this in- 
vestigation have provided a curve for CK in the range of Reynolds numbers of 
_ interest for antenna members, This curve represents the basic tool with which . 

the problem of wind-induced vibrations may be solved. In the subcritical range 
of R-values the concept of the “standard” Cx as a basis for response compu- — 
tations not only accounts for the randomness inthe amplitude of Cr, but it also 
provides independence from the mode shape of a particular vibration and the 
various other factors that affect the degree of self-amplification, = ~~ 
In the supercritical range of R-values the stochastic response concepts out- — 
lined herein identify this methodof approach tothe problem of vibrating stacks, _ 
=a and so forth, as a valid and promising means of coping with ‘past anc and 7 
_ The curves for ILF versus ASF from Phases B and C, when used in con- 
junction with the expressions for IDF, permit the computation of the degree of ; 
‘lt -amplification for any | of the types of members considered in this project. — 
Moreover, the method can be extended to other types of members by writing 
_ appropriate expressions for ASF and IDF for the particular mode shapes in- 
_ Experimental; plots of Strouhal numbers check well with previously pu published | 
iudermation, The accuracy of this agreement (within 5%) is sufficient for en- 
gineering purposes, Considering all of the possibilities for experimental er- — 
ror involved in the laboratory tests, this discrepancy is reasonable, : 

rie Effectiveness of Spoilers.—Spoiler tests show that a configuration of four 
_ helical windings is most effective in suppressing lift forces. A A spoiler diam- 

» eter of d/16 tod/8 is necessary for good suppression, and a diameter of 3 a/32- " 

reduces the lift forcesto a minimum, Although the pitchof windings is an im-_ 
portant property of the spoilers, the phenomenon is not particularly sensitive 
to the pitch. The most effective pitch in the wide range of choice (8 dto 16 d) 
was determined tobe 12d, 
= Considering the results of of spoiler optimization, it is interesting to compare ~ 
the tubular windings with a configuration of spiral fins, or “strakes,” tested by 
_C, Scruton and D, Walshe (16), They introduced a design consisting of three 
_ helical sStrakes of height up to “as and set at a pitch of 15 ad. He states that 
the strakes must be of sufficient height and number and that their helical for- 
_ mation is an essential feature. Scruton and Walshe further stated that their | q 
- effectiveness increases with size but that they need be no longer than 1/8 d. 
7 In addition, they report that the result is not sensitive to pitchand that a single . 
= strake is not effective because of the opposite- side characteristic of vortices. 
_ show characteristics similar to these strakes, with some reservations, aie 
_ Further Research.—There are two experiments of a different nature from 
this which would be of to engineers: a 


A series of with models may be conducted in a a water 
channel in order to obtain ; a continuous curve for Cx in the —/ of Reynolds 
numbers below 104, ~ 


ag = large amplitudes of oscillation with spoilers attached, whereas, vibra- aa 
tions in other types of members were readily 
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- Studies similar to those by Fung (13) and Fujino (14) should be carried 
out using large, light- -weight spring- -supported cylinders in order to 


in Lexington, Mass. 


APPENDIX L-NOTATION 


ie The following symbols, adopted for use in the paper and for the eutinane “ 
discussers, conform with "American Standards Letter Symbols for 


ican with Society representation, and approved by the 
= -cross- -sectional area, in square inches; 
Ap projected area, in square feet; 
= aeroelastic suppression factor; 


= von KarmAn coefficient, or the maximum value of the e fluctuating ‘lift: 


: oe = standard value of CK or the RMS value of CK or the ee de- 


= value of c ora barec inder; 


of c -for a with spoilers; 
constant for a a 

diameter cy linder or member, in inches; 
= modulus of elasticity | of material, in ‘pounds per square i inch; 


von Karman force = Ce q 1 Ap in pounds; 


width between the centers of two Tows of eddies behind a nonoscil-— 


lating circular cylinder, in inches, d 


= moment of inertia, ini inches; 


me incremental deflection factor; 


ring constant of an elastic system, in per 


bE 
by the M.I.T. 
| 
— 
q 
oe, 
: Fe part. 
ycles per sec- 
| 
; 
ILF 


constant a spring-supported with respect toa uni< 


form loading = the total wl, that a deflection of unity, 
pounds per inch; 
spring constant of with respect ‘to a uniform loading = the 


_ total load, w 1, that produces a maximum static deflection of unity, 

per unit length o of a spring-supported cylinder ora flexural 


member, in pound-seconds* per square inch; = 


= = stagnation pressure, , in pounds per square foots 


number = v d, in r at standard temperatur rea 


= radius o of gyration, in inc inches; 7 = 


number = (fy d d)/(17. 6 v), in air otandard temperature 


a 
spoiler ‘effectiveness; 


time, in seconds; iq 


velocity, in per hour; 


3 = critical wind velocity, at which fy = fp, in miles per hour; Ss. 
to ‘peak of vibration, in miles p per 
effective lateral loading per unit length due to vortex shedding on a 
stationary cylinder with v = = Vp» in ‘pounds: per inch; — ——— 
—_ of vibration of an elastic system, in inches; 


= of vibration at point of maximum deflection ona 
member when v* = * including of self-amplification, in 
‘ logarithmic damping decrement, or rate of decay; 


ratio of forcing function natural ‘al frequency; and 


mass density y of air, 
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THE NEUTRAL AXIS IN PLASTIC BEAMS* 


Chmene by Aris Phillips 


= 

PHILLIPS, ! ASCE. author w wishes to Abramson 
Sprague for their interest in his paper. 
a Concerning Mr. Sprague’s comments, the author believes that the applica- 
_ bility of the method presented in his paper toa beam of a given material 

depends on whether the stress-strain law of this material can be approximated 
_ with sufficient accuracy by the stress-strain law that forms the basis of the 

Tae The method has not been devised with a specific material in mind. 

It has been devised with a specific stress-strain diagram in mind, namely 
the stress- ‘Strain diagram of a perfectly plastic material. It follows that 
whether or not the ‘method is applicable to plain concrete beams depends on 
whether the perfectly plastic stress-strain diagram can be considered as a 
sufficiently good approximation of the stress-strain diagram of plain concrete. ; 
Mr. Abramson’s comments on the rotation of the neutral axis for unsym~_— 
metrical bending of beams are of much interest because they show that, even 
when the plane of bending remains constant, the neutral axis in unsymmetrical © 
bending will “not only translate but also rotate. This rotation occurs even in 
the case of a first approximation theory in which the deflections of the beam 
are assumed to have no influence on the stress distribution. In general, the 
rotation of the neutral axis should be more pronounced when the influence of 
the deflections ‘on the stress distribution is taken into account. The author 
believes that additional study of this problem, in extension of the analysis and 7 
experiments presented by Mr. , Abramson in the — sense by hi him, will 


be of considerable interest. 


October 1959, by Aris Phillips 2198), 
1 Prof., of Civ, Engrg. ‘Yale Un: Univ., » New ‘Haven, Com, 
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T. DAVISSON,® A. M. ASCE. —At least one other analog» 
‘mie of analyzing a beam on an elastic foundation is available in the litera- 
ture and is mentioned here for the sake of completeness. A. A, Wells? used 
a piano wire as the flexural member for his model studies. Although not. as. 7 
‘ae as the model presented by the author, it has the advantage of sim- 


OA review of the author’s comments concerning the solution of Eq. 1 and its — ;: 

application toa laterally loaded pile appears to be in order. By extending the - 

work of M. Hetenyi, 10 the writer!1 has presenteda derivation of the governing 
differential equation for a pile subjected to lateral, moment, and axial (w) 
ra, Eqs. 14 is the governing equation with the subscript x denoting that 

7 I, W, and Es may be functions of x. The writer used the same assumptions — - 


as the author in deriving the a ose 
an wer 
7 Eqs. 14 is nonlinear if E Sx is also made a function of y (E s y). E By solving 
the linear equation for several variations of E Sx, it is poasthle to converge 
“on a single variation that will represent E Sxy (secant modulus) for a given 
set of conditions. This is the same procedure used by the author. The axial - 
_ load term has been included for the obvious reason that a pile usually sup- 
: ports such a load as its primary function, In many cases the axial load may ; 
be ignored. However, where the soilnear the ground surface is weak or a free _ 
_— Tength of pile exists, serious errors in the flexural calculations may arise if 
axial load is ignored. . The axial load has the of magnifying moments, 
The > advantages and disadvantages of the eeveral methods of solving the 
- mens value problems represented by Eqs. 14 areof interest. Series solu- 


June 1960, by Robert Thomas (Proc, Paper 2505), 

8 Asst. Prof,, Dept, of Civ. Engrg., Univ. of [linois, Urbana, 
9 “Beams on Elastic Foundations,” ty A. A. Wells, Proceedings, Inst. of Mech, 
Engrs., Vol. 163, 1950, pp. 307-310. 
“Beams on Elastic by M, Hetenyi, University of Michigan Press, 
Ann Arbor, Michigan, 1946, 

— il “The Behavior of Flexible Verticle Piles Subjected to Moment, Shear, and Axial | 
Pere by M. T, Davisson, Thesis presented to the Univ. of Ilinois, at Urbana, Tl. in 


7 1960, in a fulfillment of the requirements for the wares of Doctor of Philosophy. 
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are wary easily obtained of Wx, and E Sx are 

Py introduced but are cumbersome and difficult to program for a digital computer. 7 
Finite difference solutions have been obtained, both by hand and digital com- — 
puter procedures, where the axial load is ignored or assumed zero. When I 
is a function of » x, the finite difference equations become quite cumbersome; a 
the situation is made worse by the addition of axial load terms. Analogs, such — 
as the mechanical methods of Wells and the author, can be quite useful. How- | 
ever, when an axial load is added, the writer envisions additional mechanical | 
errors caused by restraints at the points of load application. Boundary con- y . 

_ ditions other than free ends may also introduce mechanical difficulties. The 


writer has found that an electronic differential analyzer (analog computer) is 


readily adoptable to solving Eqs. 14, In an analog computer, the equationis 
: electrically reproduced and the accuracy of the solutions is more than satis- 


factory. Versatility in the of I W, and E with x is an n important 
| ane In summary, the writer would like to. again emphasize ' the importance | of 
7 considering the axial load in beam-columns on elastic foundations. Further- 


more, in contrast to the authors conclusions, the writer believes the analog 
_ computer isa simple and versatile means of solving Eqs.14.  — © 
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"UNDERGROUND > STRUCTURES SUBJECT AIR 


Discussion by Paul I. Rongved, and G. S. Kovacs and wis T. Frankian = 
PAUL RONGVED. 0_ The aut! authors k have ; attempted to develop an 
procedure for determining t the damage to underground structures by pseudo- 
‘steady state air overpressure. e. Unfortunately, the assumptions made for the 
behavior of the soil are, in the opinionof the writer, not realistic, and conse- 
quently, the results arrived at have little if any validity. 
7 clarify the situation the following observations are made: 


The structural ‘behavior of an underground structure is a oumgem, inter- 


5 _ Engineers have ‘not been able to ‘produce o one » general analytic procedure 
that would solve all types of structural | syatems like slabs, beams, arches, 

"shells, suspension structures,andsoon, = 
To attempt to do something - like that for a buried structure is even more. 

- difficult, but in essence this is what the author has attempted. This is believed _ 
_ to be entirely too ambitious; it should at least be deferred until well accepted 
theories for various specific buried structures have been developed. 

‘ Let us assume a uniform air overpressure varying with time applied : 

a horizontal ground surface of semi- infinite, elastic, homogenous soil. The 
vertical soil pressure under these "conditions will increase with the amount 
of the air overpressure uniformly at any point in the soil. Because we have 
a uniformly applied overpressure over uniform soil here, and the sum of all 
_ vertical forces through all horizontal sections must be ‘equal to zero at “| 
times, the above statement becomes obvious. This merely means that if ea 

- apply a uniform overpressure, such as 100 pai : at the surface, we will have a 


or any other depth below the surface. 
ws If we now assume that we have a buried structure in this soil, and the 
q structure and its foundations have the same stiffness as the soil it replaced, | 

‘the s structure will receive the same amount of pressure | as the soil had previ- 


a ously. If the structure or its foundations are harder or softer than the soil | 


they replaced, th the structure will ‘receive, correspondingly, more or less 
From these observations it is clear that our structure and its foundations 
_ should be softer than its surrounding soil. We may use a stiff structure with is 
- soft foundations or soft structure with stiff foundations or a soft structure 7 | 
with soft foundations. By the terms “soft” and “stiff” are meant a structure 
or foundation that give us a smaller or k- value 
the k-value of the surrounding soil. 
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To! ferther slain the amount of the total load carried by the ac 
7 ing soil in case of a soft rchoeieatadan , foundation on or b both, the e following observa - 
- “Most types of soil readily transmit compres 
‘siderably less shear forces, and practically no tension forces, Because ot 
this, a soil mass cannot act like a | beam or slab spanning a buried structure. — : 
The soil above a buried structure may carry its load in an infinite number ot 


= 


ways as long as it is in elastic equilibrium. Only after all possible ways of 
equilforiu the load have been exhausted will the soil reach a complete plastic 


‘Plastic equilibrium may be reached at earlier stages at individual points 
soil without causing the soil to reacha plastic 
If we can find any one aia system for the ‘soil capable of carrying | 
the load, - still allowing it to remain in its elastic equilibrium, we know that | 
its ultimate carrying capacity will be at least as great as the one we have 


found a , therefore , important to find oe beat suitable structural system 


fo 
already by ‘Engesser about 80 y yr ago, he 
_ about ‘the arching effect of soil. Arching applies to soil over linear ' 
shaped structures, For soil over circular, rectangular, and similarly 


tems in the soil that have the greatest carrying capacity determine the — 
shell structure in its” membrane condition has mainly compression 
_ forces. Soil can carry any amoun* of compression forces | up to an equivalent 
of about 1/K times the minimum side pressure exerted on the soil, K being . 
_ the coefficient of earth pressure, This holds true for soil columns, ‘soil 
arches, and soil shell structures as well. 
Looking at the example used in the article under discussion, there isa 
cubicle Structure 10 ft- «by | 10 ft-by 10 ft, buried so that the roof is 15 ft under- 
ground. No information was given on the stiffness of the foundations, but be- 
"Sse the roof will have a deflection of 1 ft at its ultimate load, it can nme 


"various structural ‘types and In the discussed article it was always: 
assumed to be a stiff cube acting as a thick slab. One can, instead, , choose 
concentrical, spherical shells that will be closer to the ‘optimum structural - 
etm. tes ten shells, one on top of the other, approximately 1 ft thick 

each, can be arranged in the soil. The shells have a radius of about q ft to 9 ft. : 

_ Following through on the analytical computation for this chosen structural 
S one sees that the carrying capacity of the soil is about at least three 
_ times greater than the value arrived at in the article discussed. This ratio 
of about 3- to-1, however, varies widely. In some ‘cases the value may be 


‘The preceding computations ‘can be made either by using a known simpli-_ 

fied, one degree of freedom , coupled dynamic design method, or by using a 
dynamic analysis v with ‘numerical integration with or without analog and digi- 

tal computers. . AS a third alternative, an analytical computation could be _ 


using simple st static- ~load design method, multiplied with an — 
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factor similar to that used for impact loads on vine cranes, and so forth, 


The impact factor must, | of course, », be developed with ‘dynamic methods orc 

The reason for the increased bearing capacity of the soil when the shell 
‘method is used can easily be understood when it is kept in mind that the shell 

action in the soilcreates horizontal components far in excess of those assumed > 
toexist inthe paper under discussion, 
7 _ The shell theory outlined has the further great advantage that the carrying _ 


"Stiffness of the structure 
3. Stiffness of the foundations of the s' 
4, _ Depth, compaction, type, and geological hist history of of the soil 


Asa a ‘result o of this it can be shown th that the dome and shell shaped under-_ 
4 ground structures induce far more carrying capacity in the overlaying soil — 


‘It is impossible to know the exact sizes and types of weapons that ca can be 


used, the” characteristics vera duration the blast wave might, therefore, 


such a way ‘that they can withstand blast waves of various duration and dis- 
tribution, This condition leads to generally stiffer structures with high buck-_ 
| ling safety; however, to insure maximum load carrying participation of the 
soil above, the foundations should be designed for the proper spring constant. 
The spring constant for soils is not exactly known, but from laboratory ex-_ 
periments and field tests values may be selected that give reasonable results, . 


Ss. KOVACS, 11 A} M. 1. ASCE and R. T. FRANKIAN, 12 M. AS ASCE.- —This 


pa 
st could be enlarged | upon by the authors in several aspects. ee 
The impression is given that the structural element is assumed to have — 
failed and then the forces that were developed in the soil are evaluated. The — 
- ‘maximum pressures are developed for failure witha rigid structure. However, | 
oe yielding necessary to reduce the applied pressures to their minimum val- 


Prides overpressure deter mined in the analysis m may be greater than n neces- 
sary for a proper analysis. 
it should also be pointed out that there is a limiting depth of burial for any 
one structure, such that additional burial will not cause a further lee 
in the overpressure felt by the structure. This limit is due to the formation 
es a ground ; arch, the height of whichis a function of the size of ‘the str structure, 
In some cases, the height of the naieal can be less than two times the width of 
asquarestructure, | 
_ ‘Further justification would be desirable for the value of 10° given as the. 
angle of internal friction for moist clays in Table 1. It appears that this was | 
only necessary so that Eq. 21 would have significance. However, the use of 
= angle of interal friction for a saturated or moist clay can be very — 


1 Intermediate Engr., Bechtel Corp., Los . Angeles, Calif, 
Project Engr., Bechtel Corp., Los Angeles, Calif. © Pa: 
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ous for cases in which the clay would be loaded rapidly and end tecntthelont time 
allowed for the excess pore-water pressure to be dissipated. tbc: 

_ The meaning of the term used for ultimate resistance, Pm, is not entirely 
clear. It was: presumed by the writers that this was the resistance to all a 
loads (static plus dynamic), However, in the example problem the ultimate. 
resistance was found insufficient to resist the static soil pressure prior to 

_ application of any dynamic load. Therefore, it must be assumed that the ulti- 
mate is that withstand dynamic load alone 


— 
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ARCH) DAM AN ‘ANA LYSIS ‘WITH A AN \N ELECTRIC 

Discussion by and Maurice L. Dickinson 


complex structures as arch dams. The implication of this method of analysis, 
so far as the domain of the civil engineer is concerned, cannot escape atten- 
tion as well as cogitation. Fortunately, the complicated as well as tedious 
analytical trial load method developed by civil | engineers, with which the 


electrical analogy may successfully compete on an economic basis, can be 


less 2 assured of success. One outstanding advantage must be conceded to the 
electrical analogy, namely the ease at which the effect of changes in assump-_ 
_ tions as to loading by water, ice, earthquake, and temperature, or as to mod- 
uli of elasticity of dam and foundations, canbe ascertained once the electrical 

_ “hookup” of the model has been cx completed, 
__ Early in 1959 the writer was called upon by Leeds, Hill, and Jewett, Con-— 

_ sulting Engineers in Los Angeles, California to participate in ‘deliberations — 
on the electrical analysis of the Blue Ridgearch dam in Arizona, that was oe : 


greatly simplified by graphic statics so that respective competition would “a 1 2 


- seribed by the author. The wealth of information on some thirty-six different _ 
“cases,” that is to say different assumptions and modifications, produced in 
- during the ten days the analog computer was occupied, was astounding. To — 

accomplish the equivalent by the analytical trial-load method would require 
many times that length of time and would certainly exceed the actual cost of 
the electrical method ¢ of analysis that, it is understood, amounted to some 
$25, 000. Of course, in civil engineering p practice the ‘design of a dam of this 
size would not be subjected to any such elaborate variations in assumptions © 
of loading, elastic properties, and structural details. Indeed, , it is not so a 
long ago that much higher arch dams were designed by the “cylinder formula” | 
based on a stress of perhaps 300 psi in a rigid arch barrel. Although these 
_ dams" still stand, such | design procedure must now be considered obsolete. 
_ ‘The author refers to tests on scale models of arch dams in Portugal that © 
“have gained due recognition by the engineering profession. Through his con-_ 
nection with Taiwan Power Company, asa member of its board of consultants, 
the: writer has recently been apprised of the results of model tests conducted | 
by the Societe’ d’ Equipement Industriel et de Laboratoirs (SEIL) in Asnierés, 
7 France on a high (760 ft +) arch dam, known as Tachien Dam, Its design was 
prepared under contract by "the engineering firm of Coyne a and Bellier in 
Paris, France, who in turn arranged the model testing with SEIL. In studying | 
- prSenEn MAES, it becomes at once ae that not only does it resemble 


“a @ August 1960, by Richard H. Mac Neal (Pree. paper 2578), 


Pox “ie ~~ ‘ 
A 4 PAUL BAUMANN,20 F, AS .— he author has presented, in a simple and _ 

a 


an Comparison Betw een Scale Model and El 


used by SEIL were as fo: follows: lows: 


(Prototype) 


Scale of model, el, S 


of loads, (Mercury) 


_ Scale of deflections, V = (16) 


— in which Em is the a Ss Modulus of the model and E,, refers to the Young’ Ss 


paris, diatomite, and water so as to closely approxi- 
_ mate Poisson’s ratio of concrete and to have a modulus of elasticity that | 
made it possible to measure strains and deformations with ease. For this 
5 purpose, 16 three-directional strain gages and deflection measuring points 
were installed on both the upstream and downstream faces. The model foun- 
— dations were cast of the same material as the model dam. This differs from 
usual prototype conditions - in that rock in place is likely to be more flexible 
than concrete. The respective discrepancy is, however, on the side of safety 
as the crown ‘stresses in an arch dam increase with Lad a in nity 
of the bed rock relative to that of concrete. 
_ The model was initially over- dimensioned relative to. the Secietitie 
et upstream and downstream faces were gradually trimmed down by means 
of a specially installed cutting machine until the final, ‘most satisfactory 
—_ was established. A certain flexibility in assumptions, although not to 
the same extent as that pertaining to the electrical model, was thereby attained, 
Contrary to the electrical analysis, no attempt was made, in the French 
> ‘model tests , to include the effect of gravity (that is weight) nor that of tem-_ 
_ perature changes. Hence, these effects had to be superposed on the “weightless” 
7 model of uniform room temperature by subsequent computations ae 
‘The reason given for “omitting a scale of weights that could easily have 
7 pom included in the computations based on structural similitude, namely be- 
cause of the uncertainty as to the effect of construction procedure on the a 
F vertical weight distribution, is not too convincing tothe writer. This, because 
the same uncertainty would seem to affect subsequent computations. The > _—, 
writer has reasons to believe that during construction a an arch dam the 


and downward, and that magnitude and direction of these forces remain | sub- 
stantially unchanged after grouting of the joints on completion of the dam. 
— the electrical analysis the effects due to gravity and temperature were 
directly ascertained by ingenious simulation as modestly | described by the 
author. On the other hand the scale model was tested as a complete shell, 
_ strictly corresponding to the prototype, rather than a combination of straight, 


rigid, horizontal arch chords and vertical beams “connected by bending-— 


— that used at the Stevenson Creek dam, described by the author, but also that 7 Ls 
_it has in prince ouch in common with the electrical “model” 
ctrical Analysis.—The scale 
= — 
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DISCUSSION 
. resisting hinges. “This articulation of an arch dam will no doubt lead to | in- a 
- creasing error in the analysis of its elastic behavior as the number of hinges — 
and intermittent chords and beams is decreased. As an absurd extreme one 


could visualize a single hinge at the crown section with rigid chords or struts 

between itandthe abutments, 

The writer feels that the distribution of load points between rigid chords 

: and beams of an arched shellis 3 a potential weakness of the electrical analysis 
that would be more or less pronounced depending on the number of load — 

_ points and intermittent members. A comparison between tests ona occas 
model and an electrical 1 ‘model, each having the same number of load points 4 

and strain gages, respectively, and in precisely the same location, would be | a 

highly recommendable because it might furnish the answer t ons that 


"MAURICE L. DICKINSON, F. ASCE.- —The significant meth- 
“ for analyzing arch dams, presented in this paper is particularly useful for 
dams located in canyons with markedly unsymmetrical abutment outlines. For 
such dams analysis, by the trial load method programmed for a digital com- 


may require simplifying ‘assumptions that could effect the validity of the re-— 


sults to an unknown, and possibly serious, degree. = OS 
_ As MacNeal points out, the operation of the computer is based on a mathe-— 


matical analogy between the behavior of the electrical components } of the 
model, and the behavior of an elastic framework representing the dam. Snel ' 
latter framework is sufficiently good representation of the dam for most de- 
sign purposes. Although there is a mathematical correspondence between the 
electrical model and the elastic framework, the method described does _ 
provide a solution of the generalized elastic equations. Thus, the method does 

é not represent an end- point: in the search for a complete solution to the problem 7 


cofarchdamanalysis, 
One of _ the principal advantages of the electric analog is the fact that it that it 


essentially represents a reasonably valid elastic model of the dam. Once the 7 yy 
electrical model has been set up, it permits easy step-by-step additions of a 


the various forces and restraints, as wellas rapid experimentation with vary-_ 


-a convenient means of investigating and demonstrating the mechanics of arch 
dam behavior i ina manner that is 3 quite useful in evaluating the various effects 

_ On the other hand, in the practical use of the electric analog, the need to _ 
“keep a considerable array of electrical components and equipment tied up | 
until the analysis is complete requires careful planning of the we! It is not 


ing boundary conditions. Consequently, the electric analog computer affords 


a. of the dectric analog. Thus, ‘whereas as the 
= ‘thickness of a particular arch ring, can be accommodated readily, changes 
in geometry of the structure must be planned in advance to avoid excessive : 
The facility with which various conditions of boundary restraint « can be a 
investigated by the electric analog leads to some important practical questions 7 


as | to how such ich degrees ¢ of restraint can be built into an actual structure. The 


{ 
ae uter can become too complex for convenient use or practicable programmin f a 
q 
— : 
pe 
Chf, Hydr, Engr., Bechtel Corp., San Francisco, Calif. 
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design and of hinged ont sliding joints in. in arch. dam structures 
would be an interesting topic for further investigation. 
Another subject that seems particularly adaptable to investigation by the | 

electric analog computer is the secondary and and more detailed analysis of 

special ‘segregated portions of an arch structure, re, such as depressed spillway 
bays, large sluice openings, -and the junction of the arch to a gravity thrust 7 
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RESTRAINED ‘COLUMN 


DONALD MANSELL. 14_The author’ pa paper is a logical of 
“the published work to which he makes reference, and in common with those - 


authors he places: severe limitations on ‘the type of problem described. The 


writer believes” that the waves f(x) can be used to solve a wider 


™ proof of the validity of the method used in deriving these 7 waves, , and in ie 

-— them, is given here, ee the general applicability of the results 


an 


(M) (M) = f; (Py) = 


October 1960, by Morris, 3, Ojalvo’ (Proc. paper ‘paper 2615). 
14 Research Student, Univ, College of Swansea, Great Brition. Bras ——o 
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i Thus, the deflected ‘curve is represented by an exact differential equation 


if the thrust-line is used as abscissa, despite the presence of sway displace-_ 
ment. In other words the solution y = f(x) is independent of sway, that may be. 
introduced into the simply as ar rotation of axeS, 


‘ 
tt is ponents to show that the ge general solution to Eq. 15 isa symmetric 
‘periodic wave if f(y) is single-valued, odd, and has positive slope. Providing 
unloading of yielded fibres does not occur, this condition is satisfactorily met 


by compression members composed of mild steel and some other structural 


materials. It is not necessary to solve the equation to demonstrate the form 
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p=. tly) dy, (16) 


3 
f(y) is od is odd, and has positive 2 ay ay is always 
_ For any value of y, for r example, a 


pm ax = - pm 


= 


always has two real vente. ‘Sinn: p_ has two zero points. 7 


_ These conditions are ‘uniquely determined because f(y) is single- valued, 


so. limit: cycle is closed and symmetric, proving that yet f(x) i isa sym- 

‘The writer has programmed th: the ir integration of Eq. 15 for a digital computer, 
‘This included the computation of the increase in centroidal ‘Strain due to 


i 
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end couples and rotations, have been compiled for use in partially- plastic 

analysis of more complex structures; this includes triangulated trusses and 

The choice of boundary conditions is of some e interest. Two of them are 
required a “complete solution of Eq. 15. The first, defining the particular 


* aq = DISCUS 
= 
Integrating Eq. 15 once, we have 
pd 
| 
(17) 
osi 
(18) 4 

a) 
cst phase diagram p v. y is symmetrical about both p and y axes. 3 
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of these functions at : any known point, as as ‘mentioned w the author. Te second 


ary condition for reasons of convenience. However, a more cogent reason for 


this choice ina general solution . may be seen from consideration of the inte- 
@ration process, 
‘The elastic solution will serve as an illustration, although the conclusion — a 
* } equally valid in the non-linear problem. The precise solution is aay 
y= -Acoskx+Bsinkx.. 
Choice of (& i 
ve of s boundary ‘condition implies that ‘A= 0. The numer 
1 integration inexact, introduces small elements of A into 


the solution. However, because the cosine term decreases in — first quadrant 


error is minimized. 


‘But Tmax is chosen, the implication is that B = 0. error 


the writer’s experience with a digital integration from ymax 
. towards y = 0 gave rise to errors in wavelength of the order of 3%, whereas, 
for- the same step-length, integration in in the Opposite direction reduced the 
error to a very small fraction of 1%. 


as 


Oa member ¢ e family of waves mav be the amplitude of either y or or one_ x 
> 
¥ 
as rounding Crrors im a COlllputer) tne roduces eleme Ol a sine 
=  termthat increases throughout the integration whereas the correct cosine term | ae 
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BEARING CAPACITY | OF ‘FLOATING ICE SHEETS* 
BEARING Cé 


(EMD, October, 1960) 

G. HOPKINS.’—Some years ago the writer was concerned with problems 

of the stress analysis of aircraft runways, and hence was familiar with the 

work of the late H. M. Westergaard. During the period from 1952 to 1954, 4 | 

_ William Prager, F. ASCE and the wetter as well as others, collaborated in ; 

the development of theories of plastic - deformations in plates under a. 


verse loads. The present approach, made on me basis of corns theory an 


a “engineering of the ‘situations involved, he is not in 
_ position to comment extensively on the nee. A particularly important ques- : 


the mechanical behavior of ice by that pay a rigid, perfectly- plastic material. 
_ Under certain conditions, such a representation of mechanical behavior was 
known to be a useful one for mild steel, a highly ductile material. . The same — 
kind of approach had been proposed by K. W. Johansen for reinforced con- 
_ crete; and now the author has proposed this approach for ice. In the develop~ — 
ment of the solutions of specific problems, it was of course necessary to 
choose a definite yield criterion, ‘and here, for reasons of mathematical sim- 
- plicity, that of Tresca had been taken, It would would be valuable to undertake 
_ experimental studies of the yield strength of sheets of ice under combined 
loading | conditions to determine the yield function as expressed in terms of 
_ principal bending moments. If it was shown to be sufficiently realistic, there 
could be theoretical advantages in working with Johansen’s (‘square’) yield 
; criterion in place of Tresca’s. 4 In addition to laboratory experimental work, 
it was necessary to undertake field experimental work to provide ao 
_ substantiation of any new theoretical approaches to such complex — as 
those with which the author was here concerned. 


a October, 1960 by G. .G. Meyerhof (Proc. paper 2627), 
Armament Research and Development Establishment, Ft. Halstead » Sevenoaks, = 
“Dynamic plastic deformations of simply-supported square plates,” 
L, L. W. Moriand, , Solids, 7, 1959, pp. 22 
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